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ABSTRACT 

Wc measure the evolution of the stellar mass function (SMF) from z — — 1 using multi- wavelength 
imaging and spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS) and the Sloan 
Digital Sky Survey (SDSS). From PRIMUS we construct an i < 23 flux-limited sample of ~ 40,000 
galaxies at z = 0.2 — 1.0 over five fields totaling sa 5.5 deg 2 , and from the SDSS we select ~ 170,000 
galaxies at z = 0.01 — 0.2 that we analyze consistently with respect to PRIMUS to minimize systematic 
errors in our evolutionary measurements. We find that the SMF of all galaxies evolves relatively little 
since z = 1, although we do find evidence for mass assembly downsizing; we measure a ~ 30% increase 
in the number density of ~ 10 10 A4q galaxies since z w 0.6, and a < 10% change in the number density 
of all > 10 11 M.q galaxies since z ~ 1. Dividing the sample into star-forming and quiescent using an 
evolving cut in specific star- formation rate, we find that the number density of ~ 10 10 M.q star-forming 
galaxies stays relatively constant since z « 0.6, whereas the space-density of > 10 11 Mq star-forming 
galaxies decreases by sa 50% between zwl and z ss 0. Meanwhile, the number density of ~ 10 10 A4q 
quiescent galaxies increases steeply towards low redshift, by a factor of ~ 2 — 3 since z sa 0.6, while 
the number of massive quiescent galaxies remains approximately constant since z~l These results 
suggest that the rate at which star-forming galaxies are quenched increases with decreasing stellar 
mass, but that the bulk of the stellar mass buildup within the quiescent population occurs around 
~ 10 10 ' 8 M.Q- In addition, we conclude that mergers do not appear to be a dominant channel for the 
stellar mass buildup of galaxies at z < 1, even among massive (> 10 11 M®) quiescent galaxies. 
Subject headings: Surveys - galaxies: evolution - galaxies: high-rcdshift - cosmology: large-scale 
structure of universe 



1. INTRODUCTION 

Surveys of the galaxy population in the nearby Uni- 
verse such as the Sloan Digital Sky Survey (SDSS; 
York et al. 2000) have found that the distribution of 
many galaxy properties, including color, morphology, 
and star formation rate (SFR), are bimodal, reflect- 
ing the existence of two broad types of galaxies: 
blue, star-forming disk galaxies, and red, quiescent 
(i.e., non star- forming) spheroidal or elliptical galax- 
ies (e.g., Kauffmann et al. 2003b; Blanton et al. 2003b; 
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Baldry et al. 2004; Wyder et al. 2007). At low red- 
shift, quiescent galaxies tend to be luminous and mas- 
sive, and are prevalent in dense environments such as 
groups and clusters, whereas star-forming galaxies typi- 
cally have lower stellar masses and are more commonly 
found in the field (Blanton & Moustakas 2009, and refer- 
ences therein). These broad empirical trends have been 
shown to persist at least to z ~ 3 (e.g., Bell et al. 2004; 
Cooper et al. 2007; Cassata et al. 2008; Brammcr et al. 
2009; Whitaker et al. 2011). Understanding how these 
two populations come into existence and evolve with cos- 
mic time, therefore, is a fundamental outstanding prob- 
lem in observational cosmology. 

Galaxy bimodality is likely a consequence of star for- 
mation in some galaxies being quenched — shut off — 
relatively quickly. A wide variety of quenching mech- 
anisms have been proposed to match the observed dis- 
tributions of galaxy properties, including: major-merger 
induced feedback from star formation and supcrmas- 
sive black holes (Springel et al. 2005; Di Mattco et al. 
2005); processes that prevent galaxies from replenish- 
ing their cold-gas supply such as virial shock heat- 
ing (Keres et al. 2005; Dekel & Birnboim 2006), often 
in concert with so-called radio-mode feedback from an 
accreting active galactic nucleus (AGN; Croton et al. 
2006; Cattaneo et al. 2006; Gabor et al. 2011); and in- 
ternal, secular quenching due to disk or bar instabil- 
ities (Cole et al. 2000; Martig et al. 2009). Late-type 
galaxies whose host halos are accreted by larger dark- 
matter halos (i.e., by groups and clusters) are also 
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susceptible to environmental quenching, such as ram- 
pressure stripping (Gunn & Gott 1972), strangulation 
or starvation (Larson et al. 1980; Balogh et al. 2000), 
and gravitational harassment (Moore et al. 1998; sec 
Boselli & Gavazzi 2006 for a review). 

By incorporating these quenching mechanisms into 
a coherent cosmological framework, modern theoreti- 
cal models have become reasonably successful at re- 
producing a wide range of galaxy properties, includ- 
ing galaxy bimodality (e.g., Hopkins et al. 2008a,b; 
Somerville et al. 2008; Bower et al. 2012). However, 
many discrepancies between observations and predictions 
persist (e.g., Fontana et al. 2006; Fontanot et al. 2009; 
Keres et al. 2009; Dave et al. 2011a,b; Lu et al. 2012; 
Weinmann et al. 2012). The exact cause of these prob- 
lems is hard to determine because of the complexities of 
modeling the small-scale physics of star formation, feed- 
back, and black-hole accretion. Consequently, empirical 
constraints on the relative fraction of quiescent and star- 
forming galaxies as a function stellar mass, redshift, en- 
vironment, AGN content, and dark-matter halo mass are 
critical for determining galaxy quenching, and the evolu- 
tion of galaxy bimodality (e.g., Weinmann et al. 2006; 
Bell et al. 2007; Kimm et al. 2009; Peng et al. 2010; 
Wetzel et al. 2012; Woo et al. 2012; Aird et al. 2012; 
Knobel et al. 2012). 

In this paper, we focus on one key aspect of this prob- 
lem by measuring the stellar mass functions (SMFs) of 
quiescent and star- forming galaxies from z = — 1, span- 
ning the last ~ 8 billion years of cosmic time. The 
SMF measures the comoving space density of galax- 
ies of a given stellar mass, making it a powerful ob- 
servational tracer of galaxy growth by in situ star for- 
mation, mergers, and galaxy transformations due to 
star formation quenching (e.g., Drory & Alvarez 2008; 
Peng et al. 2010). Measurements of the SMF are also 
important for connecting the physics of galaxy forma- 
tion to the hierarchical assembly of dark matter halos, 
and large-scale structure (e.g., Conroy & Wechslcr 2009; 
Cattaneo et al. 2011; Wang et al. 2012; Leauthaud et al. 
2012; Behroozi et al. 2012). 

Deep extragalactic surveys in the last decade have 
begun to characterize the evolution of galaxy bimodality 
over a significant fraction of the age of the Universe. At 
the highest redshifts, z > 2, studies have shown that 
although quiescent galaxies exist, they are outnumbered 
by star- forming galaxies at all stellar masses (e.g., 
Whitaker et al. 2010; Dommguez Sanchez et al. 2011). 
This early epoch is followed by a period of rapid growth 
in the space density of massive (> 10 11 Mq) quiescent 
galaxies between z ~ 2 and z ~ 1 (Arnouts et al. 2007; 
Ilbert et al. 2010; Nicol et al. 2011; Brammer et al. 
2011; Mortlock et al. 2011). By z ~ 1, the stellar 
mass dependence of galaxy bimodality as observed 
locally is largely in place: star-forming galaxies out- 
number quiescent galaxies at the low-mass end of the 
SMF, while quiescent galaxies dominate the massive 
galaxy population (e.g., Bundy et al. 2006; Borch et al. 
2006). Subsequently, between z ~ 1 and z ~ 0, the 
transformation of star-forming galaxies into quies- 
cent, passively evolving galaxies continues, leading 
to an approximately factor of two increase in the 
integrated stellar mass density of quiescent galaxies 
(Bell et al. 2004; Blanton 2006; Faber et al. 2007). 



The bulk of this stellar mass growth appears to be 
due to a rapidly rising population of intermediate- 
mass (~ 10 10 .M©) quiescent galaxies, although the 
extent to which massive galaxies also grow through 
stellar accretion (i.e., mergers) remains controversial 
(Cimatti et al. 2006; Scarlata et al. 2007; Brown et al. 
2007; Rudnick et al. 2009; Stewart et al. 2009; 
Ilbert et al. 2010; Pozzetti et al. 2010; Robaina et al. 
2010; Maraston et al. 2012). By the current epoch, qui- 
escent galaxies vastly outnumber star-forming galaxies 
above ~ 3 x 10 10 and account for more than half 

of the total stellar mass in the local Universe (Bell et al. 
2003; Baldry et al. 2004; Driver et al. 2006). 

Despite significant progress, however, the detailed evo- 
lution of the SMFs of quiescent and star-forming galaxies 
since z ~ 1 remain relatively uncertain, leaving several 
unresolved issues. One outstanding question is whether 
massive galaxies assemble their stellar mass earlier (i.e., 
at higher redshift) relative to lower-mass galaxies — that 
is, do galaxies undergo mass assembly downsizing? 12 
Mass assembly downsizing poses a significant challenge 
for theoretical models, which predict the late-time as- 
sembly of massive galaxies (i.e., mass assembly upsiz- 
ing; e.g., De Lucia et al. 2006); however, different obser- 
vational studies have reached different conclusions (e.g., 
Perez- Gonzalez et al. 2008; Fontanot et al. 2009). 

Another open question is the role of major and minor 
mergers for the stellar mass growth of galaxies at z < 1. 
Because the merger rate is notoriously difficult to mea- 
sure directly (e.g., Lotz et al. 2011), measurements of the 
SMF as a function of redshift can place complementary 
constraints on the merger-driven growth of galaxies (e.g., 
Drory & Alvarez 2008; Pozzetti et al. 2010), in addition 
to constraining the buildup of the diffuse stellar compo- 
nent (or intercluster light) of groups and clusters (e.g., 
Murante et al. 2007; Gonzalez et al. 2007). 

Finally, it remains poorly understood why the SMF of 
star-forming galaxies evolves relatively little from z = 
0—1, despite vigorous ongoing star formation. In partic- 
ular, it is not known why the stellar mass growth by in 
situ star formation balances — almost perfectly — the stel- 
lar mass growth of the quiescent galaxy population due 
to quenching (see, e.g., Arnouts et al. 2007; Martin et al. 
2007; Peng et al. 2010). 

Answers to these and related questions have remained 
elusive at intermediate redshift because a combination 
of depth, area, and large sample are needed to charac- 
terize the shape of the SMF over a large dynamic range 
of stellar mass, while simultaneously minimizing the ef- 
fects of sample variance. 13 For example, Moster et al. 
(2011) estimate a ~ 25% uncertainty in the number den- 
sity of ~ 10 11 Mq galaxies in a Az = 0.1 wide redshift 
bin at z = 0.5 due to sample variance in the ~ 2 deg 2 
Cosmic Evolution Survey (COSMOS) field, which is 
among the largest extragalactic deep fields with high- 
quality spectroscopic and medium-band photometric red- 
shifts (Scoville et al. 2007; Ilbert et al. 2009). Although 

12 Also called downsizing in stellar mass in the extensive study of 
the various manifestations of downsizing by Fontanot et al. (2009). 

13 Sample variance refers to the variation in the number density 
of galaxies along a given line-of-sight due to large-scale clustering. 
The more commonly adopted term cosmic variance should only 
strictly be used in the context of the existence of just one (observ- 
able) Universe. 
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broadband photometric rcdshifts enable the SMF to be 
constructed over larger areas (e.g., Matsuoka & Kawara 
2010), the redshift precision typically achieved by these 
methods, cr z /(l + z) ps 1% — 5%, can significantly bias 
the inferred shape of the SMF, and its evolution with 
redshift. Finally, previous studies have frequently relied 
on published measurements of the local SMF (e.g., from 
the SDSS), although the amount of evolution inferred 
can be significantly affected by systematic differences in 
how stellar masses are derived (Marchesini et al. 2009; 
Bernardi et al. 2010). 

We alleviate many of these issues by measuring the 
evolution of the SMF at intermediate redshift using 
data from the PRism MUlti-object Survey (PRIMUS; 
Coil et al. 2011; Cool et al. 2013). From PRIMUS we 
select ~ 40, 000 galaxies at z = 0.2 — 1 with high-quality 
spectroscopic redshifts and deep multi-wavelength imag- 
ing in the ultraviolet (UV) from the Galaxy Evolution Ex- 
plorer (GALEX; Martin et al. 2005), in the mid- infrared 
from the Spitzer Space Telescope (Werner et al. 2004) In- 
frared Array Camera (IRAC; Fazio et al. 2004), and in 
the optical and near-infrared from a variety of ground- 
based surveys. The broad wavelength coverage allows us 
to estimate stellar masses and SFRs using detailed spec- 
tral energy distribution (SED) modeling, and to robustly 
identify quiescent and star-forming galaxies over the full 
redshift range. The total area subtended by this sam- 
ple is ps 5.5 deg 2 , making it the largest statistically com- 
plete sample of faint galaxies with spectroscopic redshifts 
ever assembled. Furthermore, we construct the SMF at 
z ss 0.1 using a sample of ~ 170,000 SDSS galaxies at 
z = 0.01 — 0.2 with UV, optical, and near-infrared pho- 
tometry from GALEX, the Two Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006), and the Wide-field In- 
frared Survey Explorer (WISE; Wright et al. 2010) that 
we analyze using the same methodology as PRIMUS to 
minimize systematic errors in the evolutionary trends we 
measure. 

Using these data, we measure the evolution in the num- 
ber and stellar mass density of quiescent and star-forming 
galaxies since z ps 1, in order to quantify their stellar 
mass growth by star formation and mergers, and to con- 
strain the physical mechanisms responsible for quench- 
ing. In Section 2 we present the GALEX, optical, and 
IRAC photometry of galaxies in the PRIMUS fields that 
we use, and we describe how we construct our local SDSS 
comparison sample. We select our parent sample of qui- 
escent and star-forming galaxies in Section 3, and in Sec- 
tion 4 we describe the methodology we use to construct 
a statistically complete SMF for both quiescent and star- 
forming galaxies in seven redshift bins from z = 0—1. We 
present our SMFs and quantify the number and stellar 
mass density evolution of each galaxy type in Section 5, 
and quantify the role of galaxy growth by mergers and 
star-formation quenching in Section 6. Finally, we sum- 
marize our results in Section 7. 

Given the length of the paper, readers interested in 
our principal results can skip ahead to Section 5 (see 
especially Figures 11 and 12), and to the interpretation 
of our results in Section 6. 

We adopt a concordance cosmology with Q m = 0.3, 
Qa = 0.7, and h 70 = H a /(70 km s" 1 Mpc" 1 ) = 1.0, and 
the AB magnitude system (Oke & Gunn 1983) through- 



out. Unless otherwise indicated, all stellar masses and 
SFRs assume a universal Chabrier (2003) initial mass 
function (IMF) from 0.1 - 100 M Q . 

2. OBSERVATIONS 

Our analysis of the SMF at intermediate redshift com- 
bines multi-wavelength imaging of five distinct extra- 
galactic deep fields with spectroscopic redshifts from 
PRIMUS. In Section 2.1 we briefly describe PRIMUS 
and introduce the fields that we analyze. We describe 
our analysis of the deep GALEX /UV imaging of these 
fields in Section 2.2, and the optical and mid-infrared 
photometric catalogs we use in Section 2.3. Finally, in 
Section 2.4 we describe how we construct our z ps 0.1 
SDSS galaxy sample. 

2.1. PRIMUS 

PRIMUS is a large faint-galaxy intermediate-redshift 
survey which obtained precise (cr z /(l + z ) ~ 0.5%) spec- 
troscopic redshifts for a statistically complete sample 
of ~ 120,000 galaxies at z ps - 1.2. The survey 
targeted galaxies in seven distinct extragalactic deep 
fields, totaling ~ 9 deg 2 , with a wealth of ancillary 
multi-wavelength imaging from the X-ray to the far- 
infrared. PRIMUS was conducted using the IMACS 
spectrograph (Bigelow & Dressier 2003) on the Magel- 
lan I Baade 6.5 m telescope with a slitmask and low- 
dispersion prism. This novel experimental design yielded 
low- resolution (A/AA ~ 40) spectra for ~ 2000 objects 
per slitmask, which is a factor of ~ 10 higher multi- 
plexing rate than traditional spectroscopic redshift sur- 
veys (see also Kelson et al. 2012). The full details of 
the survey design, targeting, and data summary are in 
Coil et al. (2011), while the details of the data reduction, 
redshift fitting and precision, and survey completeness 
can be found in Cool et al. (2013). 

In this paper we restrict our analysis to the fields 
targeted by PRIMUS with GALEX and Spitzer /IRAC 
imaging. Three of these fields are part of the Spitzer 
Wide-area Infrared Extragalactic Survey (SWIRE 14 ; 
Lonsdale et al. 2003): the European Large Area ISO 
Survey - South 1 field (ELAIS-S1 15 ; Oliver et al. 2000); 
the Chandra Deep Field South SWIRE field (hereafter, 
CDFS); and the XMM Large Scale Structure Survey field 
(XMM-LSS; Pierre et al. 2004). In detail, the XMM-LSS 
field in PRIMUS consists of two separate (but spatially 
adjacent) subfields: the Subaru/XMM-Newton DEEP 
Survey field (XMM-SXDS 16 ; Furusawa et al. 2008), and 
the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) field (hereafter, XMM-CFHTLS 17 ). These 
fields were targeted by PRIMUS using different photo- 
metric catalogs; therefore, because of the slightly differ- 
ent selection functions we treat them separately in our 
analysis of the SMF. Finally, we include in our analysis 
the well-studied COSMOS 18 field (Scovillc et al. 2007), 
giving a total of five distinct fields. 

2.2. Ultraviolet Photometry 

14 http:/ /swirc. ipac.caltcch.edu/swire/swire. html 

15 http:/ /dipastro. pd.astro.it/esis 

16 http:/ /www. naoj.org/Sciencc/SubaruProject/SDS 

17 http:/ /www. cfht.hawaii.edu/Scicnce/CFHLS 

18 http:/ /cosmos. astro. caltech.edu 
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All five fields in our sample were observed in the 
near-UV (NUV) and far-UV (FUV) as part of the 
GALEX Deep Imaging Survey (DIS; Martin ct al. 2005; 
Morrissey et al. 2005). The mean exposure times of these 
observations ranges from 58—78 ks in our COSMOS field, 
to between 15 - 42 ks in our ELAIS-S1, CDFS (but one 
with 90 ks), and XMM-SXDS and XMM-CFHTLS (but 
one with 150 ks) fields, making these among the deepest 
UV observations ever obtained. 

Given the 4'/2 and 5'/3 FWHM point-spread function 
(PSF) of the GALEX telescope in the FUV and NUV, 
respectively, in > 10 ks exposures source confusion is a 
significant problem, requiring great care when extracting 
a photometric catalog. Therefore, to minimize contam- 
ination from neighboring sources, we use the Bayesian 
photometric code EMphot, which is based on the expec- 
tation maximization (EM) algorithm of Guillaumc et al. 
(2006). EMphot uses optical positional priors to mea- 
sures the UV fluxes of objects in the GALEX images by 
adjusting a model of the GALEX PSF. The prior posi- 
tions are based on deep, high-resolution optical imaging 
of the same field in the bluest available band (see Sec- 
tion 2.3). EMphot has been used successfully in sev- 
eral previous studies to analyze deep GALEX imaging 
(Xu ct al. 2005; Zamojski et al. 2007; Salim ct al. 2009; 
Hammer ct al. 2010), so we refer the interested reader to 
those papers for additional details. 

2.3. Optical & Mid-Infrared Photometry 

In this section we describe the ground-based opti- 
cal and Spitzer /IRAC mid- infrared photometric catalogs 
that we utilize in each field. As most of these catalogs 
are publicly available, we defer many of the finer details 
to the corresponding data release documentation and pa- 
pers cited below; additional details can also be found in 
Coil et al. (2011). 

Our general strategy for obtaining integrated (total) 
fluxes in all photometric bands is to use circular aper- 
ture photometry to constrain the shape of the SEE) of 
each galaxy, and then to tie the overall normalization 
of the SED to an estimate of the total magnitude in 
the detection band (usually i' or R). The strengths of 
this procedure are that aperture colors are less affected 
by neighboring sources, and they typically have higher 
signal-to-noise ratios than total magnitudes measured in 
each band independently. In three of our fields — XMM- 
SXDS, XMM-CFHTLS, and COSMOS— aperture colors 
have been measured from point spread function (PSF) 
matched images and so should be very accurate; in the 
other two fields— CDFS and ELAIS-S1— the aperture 
colors are not measured from PSF-matchcd mosaics, so 
for these sources we adopt slightly higher minimum pho- 
tometric uncertainties when performing our SED model- 
ing (see Section 4.1). 

2.3.1. CDFS, ELAIS-S1, XMM-SXDS, and 
XMM-CFHTLS 

In the CDFS field we use the ground-based optical 
photometric catalogs distributed as part of the SWIRE 
Data Release 2/3, as described in the SWIRE Data De- 
livery Document. 19 The SWIRE team obtained Ug'r'i'z' 

19 http://irsa.ipac.caltech.edu/data/SPITZER/SWIRE 



imaging of this field using the MOSAIC-II imager at the 
CTIO/Blanco 4 m telescope, achieving a depth in each 
band of 25.2, 25.3, 25.2, 24.4, and 23.8 mag (5cr) for 
a point source in a 2" diameter aperture (Norris ct al. 
20 06). 20 Source catalogs were generated (by the SWIRE 
team) in each bandpass individually using the Cambridge 
Astronomical Survey Unit (CASU 21 ) pipeline, and then 
merged using the Spitzer Science Center bandmerge 22 
software package. We adopt the fluxes measured in a 
fixed 2'.'4 diameter aperture, scaled to match CASU's es- 
timate of the integrated i' flux. 

In the ELAIS-S1 field we use the BVR and Iz cat- 
alogs published by Berta et al. (20 06 , 20 08). 23 The 
BVR imaging was obtained using the Wide Field Im- 
ager (WFI) at the 2.2 m La Silla ESO-MPI telescope, 
achieving a depth of 24.9, 25.0, and 24.7 mag in B, 
V, and R, respectively. The /- and z-band observa- 
tions were obtained using the Visible Multi Object Spec- 
trograph (VIMOS; Le Fevre et al. 2003) camera at the 
VLT 8.2 m telescope, reaching a depth of ~ 23.3 and 
~ 22.8 mag, respectively. All the quoted depths corre- 
spond to the 95% completeness limit for point sources. 
The publicly released photometric catalogs were gener- 
ated using SExtractor (Bertin & Arnouts 1996). We 
adopt the aperture fluxes measured in a 3'.' 3 diameter 
aperture, scaled to the MAG_AUTO (total) i?-band mag- 
nitude. 

Finally, in the XMM-SXDS and XMM-CFHTLS fields 
we rely on the high-fidelity u*g'r'i'z' photometric cat- 
alogs generated as part of the CFHTLS Archive Re- 
search Survey (CARS; Erben et al. 20 09). 24 These cat- 
alogs are based on deep imaging obtained as part of the 
CFHTLS- Wide survey using the CFHT/Megacam cam- 
era (Boulade et al. 20 03). 25 The CARS mosaics reach a 
5cr depth of 25.2, 25.3, 24.4, 24.7, and 23.2 mag in a 2" di- 
ameter aperture in u* , g' , r , i' , and z', respectively. The 
CARS photometric catalogs were generated from PSF- 
matched images using SExtractor in dual-image mode, 
with the unconvolved i' mosaic serving as the detection 
image. We use the fluxes measured in a 3" diameter 
aperture, scaled to the MAG_AUTO magnitude measured 
from the unconvolved i'-band mosaic. 

In addition to the ground-based optical observations 
described above, all four of the preceding fields were also 
observed at 3.6, 4.5, 5.8, and 8 fim with Spitzer /IRAC 
as part of SWIRE. We use the SWIRE Data Release 
2/3 IRAC catalogs matched to the galaxies in PRIMUS 
using a 1" search radius. These catalogs are complete 
for point sources to a 5a depth of 22.2, 21.5, 19.8, and 
19.9 mag at 3.6, 4.5, 5.8, and 8 /im, respectively. Fol- 
lowing the SWIRE Data Delivery document, we use the 
fluxes measured in a 3'.' 8 diameter circular aperture, mul- 
tiplied by an aperture correction derived in each band 
from isolated point sources (see also Surace et al. 2004; 
Ilbert et al. 20 09). 26 These aperture-corrected fluxes op- 

20 http: //www. astro . caltech. edu/$\sim$bsiana/cdf s_opt 

21 http://casu.ast.cam.ac.uk 

22 http: // ssc . spitzer . caltcch.edu / dataanalysistools / tools /bandmerge 

23 http:/ /www. astro. unipd.it/esis 

24 ftp: / /marvinweb. astro. uni-bonn.de/data_products/CARS_catalogues 

25 http://www.cfht.hawaii.edu/Science/CFHLS 

26 p or re f erencej the aperture correction factors we use, as mea- 
sured by the SWIRE team, are 0.736, 0.716, 0.606, and 0.543 at 
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timize the trade-off between crowding, which favors a 
smaller aperture, and signal-to-noise ratio, which favors a 
larger aperture, and provide a reasonably accurate mea- 
surement of the total IRAC flux (see the SWIRE Data 
Delivery Document). 

2.3.2. COSMOS 

In the COSMOS field we rely on the multi-wavelength 
photometric catalog publicly released in April 2009 by 
the COSMOS team (see Capak et al. 2007). 27 This 
catalog includes Vjg + r + i + z + imaging obtained us- 
ing the Suprime-Cam instrument (Miyazaki et al. 2002) 
on the Subaru 8.2 m telescope; u*i* imaging from 
the Megacam camera (Boulade et al. 2003) on the 
3.6 m Canada- France-Hawaii Telescope (CFHT); and 
-Ks-band imaging obtained using the Wide-field InfraRcd 
Camera (WIRCam; Puget et al. 2004) on the CFHT 
(McCracken et al. 2010). 28 In a 3" diameter aperture, 
these images achieve a 5<r depth of 26.5, 26.6, 26.1, 23.5, 
and 25.1 mag in u*, g + , Vj, r + . i + , i* , and z + , re- 
spectively; the near-infrared imaging is more than 90% 
complete for point sources to K s = 23. The COSMOS 
optical/near-infrared photometric catalog was generated 
from PSF-matchcd mosaics using SExtractor in dual- 
image mode, with the unconvolved i + mosaic as the de- 
tection image. We use the fluxes measured in a 3" di- 
ameter aperture, scaled to the i + MAG_AUTO magnitude 
measured from the unconvolved z + -band mosaic. 

We supplement these optical/near- infrared data with 
mid-infrared photometry generated by Mendez et al. 
(2013) using the public Spiteer-COSMOS 29 (S- 
COSMOS) IRAC mosaics (Sanders et al. 2007). The 
catalogs were generated using the same tools utilized by 
the SWIRE team (see Section 2.3.1), thereby ensuring 
that our IRAC photometry is consistent across all five 
fields; nevertheless, Mendez et al. (2013) demonstrate 
that their photometric measurements and uncertainties 
agree in the mean with the S-COSMOS/IRAC catalogs 
publicly released by the COSMOS team in June 2007. 30 
The catalogs arc statistically complete for point sources 
brighter than 24.0, 23.3, 21.3, and 21.0 mag (5tr) in each 
of the four IRAC channels. As for the other four fields, 
we use the aperture-corrected 3'.'8 diameter aperture 
fluxes and match each IRAC source to our PRIMUS 
sample using a 1" search radius. 

2.4. SDSS-GALEX Sample 

An accurate measurement of the SMF at the current 
epoch is crucial because it provides a low-redshift an- 
chor against which we can quantify the stellar mass 
buildup of galaxies through cosmic time. Although 
many previous studies have measured the local SMF 
for the global galaxy population (see Section 5.1), few 

3.6, 4.5, 5.8, and 8 (im, respectively. (The fluxes are divided by 
these factors.) 

27 http://irsa.ipac.caltech.edu/data/COSMOS 

28 We do not use the Subaru Bj- and Vj-band imaging of 
the COSMOS field, nor the UKIRT J-band imaging, because 
of their larger-than-average photometric zeropoint uncertainties 
(Ilbert et al. 2009). 

29 http://irsa.ipac.caltech.edu/data/SPITZER/S-COSMOS 

30 http:/ /irsa. ipac.caltech.edu/data/COSMOS/gator_docs/ 
scosmos_irac_colDcscriptions.html 



have exploited the bimodality of the SFR versus stel- 
lar mass diagram to separately measure the SMFs of 
quiescent and star-forming galaxies. Moreover, given 
the susceptibility of stellar mass estimates to a vari- 
ety of model-dependent systematic uncertainties (see, 
e.g., Marchesini et al. 2009; Bernardi et al. 2010; Ap- 
pendix B), it is important that we measure the local 
SMF using the same assumptions and methodology used 
to generate the SMF of intermediate-redshift galaxies. 

With these considerations in mind, we select a sam- 
ple of low-redshift galaxies using the SDSS Data Re- 
lease 7 (DR7; Abazajian et al. 2009), which provides 
high-fidelity ugriz photometry and spectroscopic red- 
shifts for hundreds of thousands of galaxies in the nearby 
Universe. Specifically, we select 504, 437 galaxies from 
the New York University Value-Added Galaxy Catalog 
(VAGC 31 ; Blanton et al. 2005b) that satisfy the main 
sample criteria defined by Strauss et al. (2002), and have 
Galactic extinction corrected (Schlegel et al. 1998) Pet- 
rosian magnitudes 14.5 < r < 17.6, and spectroscopic 
redshifts 0.01 < z < 0.2. Excluding areas of the 
survey that were masked because of bright stars and 
other artifacts (Blanton et al. 2005b), this sample covers 
6956 dcg 2 . The VAGC also provides an estimate of the 
statistical weight for each galaxy, which we use to correct 
the sample for targeting incompleteness (generally due 
to fiber collisions) and redshift failures (Blanton et al. 
2003a). 

Next, we restrict the parent SDSS sample to galax- 
ies with medium-depth observations from GALEX; UV 
photometry is needed to effectively divide the sam- 
ple into quiescent and star-forming galaxies (see Sec- 
tion 3.2). First, we retrieve the positions of all the 
GALEX tiles publicly available as part of the GALEX 
Release 6 (GR6) with a total exposure time greater 
than 1 ks, which consists of more than 5400 tiles cov- 
ering roughly 4450 deg 2 of the sky. The 1 ks expo- 
sure time cut is the minimum time necessary to ob- 
tain an accurate measurement of the total UV flux of 
extended sources at the typical redshift of galaxies in 
the SDSS (Wyder et al. 2007). Next, we use mangle 32 
(v2.2; Hamilton & Tegmark 2004; Swanson et al. 2008) 
to construct the angular selection function of the joint 
SDSS- GALEX sample. For the SDSS we use the mangle 
polygons distributed as part of the VAGC, while for each 
GALEX pointing we adopt a simple 1?1 diameter circu- 
lar field-of-view (Morrissey et al. 2007). The final sample 
comprises 169, 727 SDSS galaxies with GALEX imaging 
distributed over 2505 deg 2 . 

With the relevant list of objects in hand, we use the 
MAST/CasJobs 33 interface and a 4" diameter search ra- 
dius (Budavari et al. 2009) to retrieve the NUV and FUV 
photometry of the galaxies in our sample. We resolve 
duplicate GALEX observations due to overlapping tiles 
by selecting the measurement with the highest signal-to- 
noise ratio, and adopt the SExtractor MAG_AUTO mag- 
nitude (Bertin & Arnouts 1996) measured in the NUV 
band as an estimate of the total NUV flux. For the FUV 
photometry we use the FUV flux measured inside the 
elliptical aperture defined in the NUV band to ensure 

31 http:/ /sdss. physics. nyu.edu/vagc 

32 http:/ /space. mit.edu/~molly/mangle 

33 http:/ /galex. stsci.edu/casjobs 
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accurate galaxy colors. 

For the ugriz bands we use the SDSS model mag- 
nitudes, which provide reliable, high signal-to-noise 
measurements of the optical colors of each galaxy 
(Stoughton et al. 2002). We scale the ugriz photome- 
try to the r-band cmodel magnitude, which provides the 
most reliable estimate of the total (integrated) galaxy 
flux irrespective of galaxy type (Abazajian et al. 2004; 
Bernardi et al. 2010; Blanton et al. 2011). Finally, we 
supplement our UV and optical photometry with inte- 
grated JHK S magnitudes from the 2MASS Extended 
Source Catalog (XSC; Jarrett et al. 2000), and with in- 
tegrated photometry (or upper limits) at 3.4 and 4.6 /im 
from the WISE All-Sky Data Release. 34 

3. SAMPLE SELECTION 

In Section 3.1 we define the PRIMUS parent sample, 
and in Section 3.2 we describe the criteria we use to select 
quiescent and actively star-forming galaxies as a function 
of rcdshift. 

3.1. Parent Sample 

We select our parent sample from the statistically com- 
plete primary sample of galaxies observed by PRIMUS 
(see Coil et al. 2011 for full details). In Table 1 we list 
the optical selection band and the corresponding mag- 
nitude limits we impose to define the parent sample in 
each of our five fields. These limits are 18 < i 1 < 23 
in the CDFS, XMM-SXDS, and XMM-CFHTLS fields, 
18 < I < 23 in COSMOS, and 18 < R < 23.2 in ELAIS- 
Sl. Due to differences in the PRIMUS experimental de- 
sign, the sample of galaxies targeted for spectroscopy in 
the CDFS and ELAIS-S1 fields were also required to be 
detected in IRAC imaging at 3.6 /jm; therefore, we fur- 
ther require our parent sample of galaxies in these two 
fields to have 17 < [3.6] < 21. We emphasize that hav- 
ing our sample flux-limited in both the optical and mid- 
infrared does not preclude us from determining the SMF 
in these two fields, as we account for both flux limits in 
our analysis (see Section 4.3). 

Next, we construct the angular selection (window) 
function of our sample by joining the PRIMUS, GALEX, 
and IRAC window functions using mangle. The final 
window function includes regions of the sky with GALEX 
and IRAC imaging, and ensures coverage from two or 
more PRIMUS slitmasks, thereby minimizing targeting 
incompleteness due to slit collisions. The solid angle 
of each field ranges from 0.80 deg 2 in ELAIS-S1, to 
1.70 deg 2 in the XMM-CFHTLS field, totaling 5.50 deg 2 
(see Table 1). 

Finally, we select objects spectroscopically classified 
as galaxies by the PRIMUS pipeline (thereby excluding 
stars and broad-line AGN; see Cool et al. 2013 for de- 
tails) with high-quality {Q > 3) spectroscopic redshifts 
in the range z = 0.2 - 1.0. Below z = 0.2 PRIMUS is 
severely limited by sample variance, while our upper red- 
shift cut eliminates < 5% of the primary sample. The 
rcdshift confidence cut Q > 3 balances the need for a 
large sample while minimizing the catastrophic outlier 
rate and maximizing the redshift precision. Based on a 
comparison with redshifts derived from high-resolution 

34 http: //wise2 . ipac . caltech . edu/docs/release/allsky 



spectroscopy, we estimate that the redshift precision of 
our sample is cr z /(l + z) w 0.4%, with a catastrophic 
outlier rate of < 3%. Table 1 lists the final number of 
galaxies in our sample. 

To correct our sample for targeting incompleteness and 
rcdshift failures, we use the statistical weights described 
by Coil et al. (2011) and Cool et al. (2013). Briefly, we 
assign a statistically derived weight, u>i, to each galaxy 
given by 

(/target X /collision X /success) 3 (-0 

where /target is the fraction of galaxies that passed the 
PRIMUS magnitude- and density-dependent target se- 
lection criteria; /collision is the fraction of potential tar- 
gets observed spectroscopically (i.e., whose spectrum 
would not collide, or overlap with another potential tar- 
get); and /success is the fraction of galaxies in a given bin 
of apparent magnitude and color (typically g— r or B—R) 
that yielded a robust (i.e., Q > 3) redshift. Both the tar- 
geting fraction (/target ~ 80%; Coil et al. 2011) and the 
fraction of targets observed (/collision ~ 95%; Cool et al. 
2013) are high in PRIMUS due to its novel experimental 
design and survey strategy. The redshift success rate de- 
creases smoothly with apparent magnitude — essentially, 
spectroscopic signal-to-noisc ratio — and varies weakly 
with observed-frame color. For reference, / SUC ccss de- 
creases from > 75% at i = 21, to ps 45% at i = 22.5, and 
to w 30% at the limit of our survey, i ss 23 (Cool et al. 
2013). We emphasize that the lack of a significant trend 
of /success with observed-frame color, as well as exten- 
sive comparisons of the PRIMUS redshifts with redshifts 
derived from high-resolution spectroscopic surveys for 
both intrinsically red and blue (i.e., quiescent and star- 
forming) galaxies, indicate that the redshift success in 
PRIMUS is not a strong function of galaxy type. 

3.2. Selecting Quiescent and Star-Forming Galaxies 

A variety of techniques have been proposed to sepa- 
rate passively evolving galaxies from galaxies with on- 
going star formation (see, e.g., Williams et al. 2009; 
Pozzetti et al. 2010), but at a basic level all methods ex- 
ploit to varying degrees the existence of galaxy bimodal- 
ity (see Section 1). 

Here, we leverage our broad wavelength coverage and 
precise spectroscopic redshifts (see Section 2) to measure 
accurate stellar masses and SFRs for the galaxies in our 
sample using iSEDfit, a new Bayesian SED modeling 
code (see Section 4.1 for details). With these quanti- 
ties in hand, we divide the galaxy population into star- 
forming and quiescent based on whether they lie on or 
below the so-called star formation (SF) sequence. The 
SF sequence (also called the main sequence of star for- 
mation; Nocskc et al. 2007) is the correlation between 
SFR and stellar mass exhibited by star-forming galaxies 
at least to z ~ 2 (e.g., Oliver et al. 2010; Karim et al. 
2011). In Figure 1 we plot SFR versus stellar mass in 
seven rcdshift bins from z = — 1 for both our SDSS- 
GALEX and PRIMUS samples. We find a well-defined 
SF sequence whose amplitude increases smoothly toward 
higher rcdshift, and a distinct population of quiescent 
galaxies that fall below the SF sequence at a given stel- 
lar mass (e.g., Salim et al. 2007; Elbaz et al. 2007). We 
postpone a more detailed discussion of the evolution of 



Evolution of the Stellar Mass Function 



7 



z = 0.20-0.30 



e 

s - 



K 
m 





z = 0.40-0.50 





z = 0.65-0.80 




OS 
fe- 
rn 




9 10 11 

log (M / M s ) 



12 



10 11 12 
log (M / M e ) 



Fig. 1. — Star formation rate (SFR) vs. stellar mass in seven bins 
of rcdshift from z = — 1 based on our SDSS- GALEX (upper-left 
panel) and PRIMUS (subsequent six panels) samples. We divide 
our sample into star-forming or quiescent according to whether 
they lie above or below the dashed line, respectively; this line is 
parallel to the star formation (SF) sequence at z 0.1 and evolves 
with rcdshift according to equation (2). 

the SF sequence to z = 1 using PRIMUS to another 
paper. 

To divide the galaxy population we use an evolving 
cut that traces the lower envelope of the SF sequence in 
each redshift bin. In detail, we first rotate the SFR versus 
stellar mass diagram using the power-law slope of the SF 
sequence derived by Salim et al. (2007); Salim et al. find 
SFR oc .A/f ' 65 for galaxies at z « 0.1, which is also a good 
fit to our SDSS- GALEX sample. Next, we construct 
the histogram distribution of "rotated" SFRs, given by 
log (SFR rot ) = log (SFR) - 0.65(log M - 10), where SFR 
is in units of A4q yr _1 and M. is in M.q. and identify 
(by eye) the minimum of the bimodality in each redshift 
bin. Finally, we fit the minimum of the bimodality versus 
redshift to obtain the following linear relation: 

log(SFR min ) = -0.49 + 0.65 log (Af- 10) + 1.07 (z-0.1). 

We classify each galaxy in our sample into star-forming 
and quiescent based on whether its SFR and stellar mass 
place it above or below the SFR given by equation (2), 
interpolated at the redshift of the galaxy. 

4. BUILDING THE STELLAR MASS FUNCTION 

In this section we describe how we infer the stellar 
masses and SFRs for the galaxies in our sample (Sec- 



tion 4.1), present the technique we use to construct a 
non-parametric estimate of the SMF (Section 4.2), and 
calculate the stellar mass above which our full (hereafter, 
the all sample), quiescent, and star- forming galaxy sam- 
ples are statistically complete as a function of redshift 
(Section 4.3). 

4.1. Stellar Masses and Star Formation Rates 

Modeling the broadband SEDs of galaxies using stel- 
lar population synthesis models has become a power- 
ful technique for inferring their physical properties (see 
the recent review by Walcher et al. 2011, and refer- 
ences therein). We have developed iSEDfit, a suite 
of routines written in the IDL programming language 
to determine within a simplified Bayesian framework 
the stellar masses, SFRs, and other physical properties 
of galaxies from their observed broadband SEDs (e.g., 
Kauffmann et al. 2003a; Salim et al. 2007; Auger et al. 
2009). We describe iSEDfit in more detail in Ap- 
pendix A, but in essence the code uses the redshift and 
observed photometry of each galaxy to compute the sta- 
tistical likelihood of a large ensemble of model SEDs — 
generated using population synthesis models — spanning 
a wide range of observed colors and physical proper- 
ties (stellar mass, age, metallicity, star formation his- 
tory, dust content, etc.). Random draws of the model 
parameters are chosen from user-defined prior parame- 
ter distributions using a Monte Carlo technique. Once 
the posterior probability distribution function (PDF) has 
been computed, the marginalized PDF of the quantity 
of interest, such as the stellar mass, follows from the 
probability-weighted histogram of the corresponding pa- 
rameter values. The median of the posterior PDF can 
then be adopted as the best estimate of that parameter, 
and the uncertainty can be derived from the cumulative 
distribution function (see Appendix A for more details). 

Although SED modeling is conceptually straightfor- 
ward, the inferred physical properties depend on which 
population synthesis models and prior parameters are 
adopted (e.g., Marchesini et al. 2009; Kajisawa et al. 
2009; Muzzin et al. 2009). For example, differences 
in the stellar libraries (e.g., theoretical vs. empirical) 
among population synthesis models, and the exact 
treatment of post-main sequence stellar evolution can 
result in widely different predictions of the time- and 
metallicity-dependent spectral evolution of even simple 
(i.e., coeval) stellar populations (Conroy et al. 2009; 
Conroy & Gunn 2010; Mancone & Gonzalez 2012). 
Moreover, the exact choice of prior parameters (e.g., 
dust attenuation curve, bursty vs. smooth star formation 
histories, treatment of metallicity evolution, etc.) can 
also have a significant effect on the derived properties 
(Kannappan & Gawiser 2007; Perez- Gonzalez et al. 
2008; Carter et al. 2009; Longhetti & Saracco 
2009; Muzzin et al. 2009; Kajisawa et al. 2009; 
Marchesini et al. 2009; Maraston et al. 2006, 2010; 
Conroy et al. 2010; Pforr et al. 2012). Finally, SED 
modeling typically includes an implicit assumption of a 
fixed, universal IMF — an IMF that does not vary with 
rcdshift or galactic physical conditions. Throughout this 
paper wc make the same simplifying assumption, al- 
though recent observations suggest the IMF may not be 
as universal as once was thought (Dave 2008; Trcu et al. 
2010; van Dokkum et al. 2010; Cappellari et al. 2012, 
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but see Bastian et al. 2010 for a critical review of the 
evidence). 

An exhaustive investigation of the preceding issues is 
beyond the scope of the current study. Nevertheless, we 
would like to have a qualitative and quantitative sense of 
which of our results are (in)sensitive to the exact choice 
of population synthesis models and priors. Therefore, we 
proceed by presenting our principal results in the main 
body of the paper using a fiducial set of SED modeling 
assumptions, and in Appendix B we examine the effect 
of varying these assumptions on our conclusions. 

Here, we briefly summarize the default population syn- 
thesis models and prior parameters we use, but refer 
the interested reader to Appendices A and B for addi- 
tional parameter definitions and details. Our fiducial 
stellar masses and SFRs are derived using the Flexi- 
ble Stellar Population Synthesis (FSPS 35 ) models (v2.3; 
Conroy et al. 2009; Conroy & Gunn 2010), based on the 
Chabrier (2003) IMF from 0.1 - 100 M e . We con- 
sider exponentially declining star formation histories 
with stochastic bursts of varying onset, strength, and du- 
ration superposed (Kauffmann et al. 2003a; Salim et al. 
2007), and allow a wide range of galaxy ages and possi- 
ble star formation histories. Finally, we assume sensibly 
distributed priors on stellar mctallicity and dust attenu- 
ation, and adopt the time-dependent attenuation curve 
of Chariot & Fall (2000). 

The photometric bands we use vary with the sample. 
In our SDSS- G ALEX sample we fit to 12 bands of pho- 
tometry: FUV and NUV from GALEX; ugriz from the 
SDSS; JHK 8 from 2MASS; and the 3.4 and 4.6 /im bands 
from WISE (see Section 2.4). In our PRIMUS sample 
we fit to our GALEX FUV and NUV photometry (see 
Section 2.2), the two shortest IRAC bands at 3.6 and 
4.5 /im, and five optical bands, except in COSMOS where 
we fit to seven optical and near-infrared bands (see Sec- 
tion 2.3). We do not fit to the two longer- wavelength 
IRAC channels at 5.8 and 8 /im because of the poten- 
tial contributions from hot dust and polycyclic aromatic 
hydrocarbon (PAH) emission lines at these wavelengths 
(e.g., Smith et al. 2007), which iSEDf it does not cur- 
rently model. When fitting, we assume minimum photo- 
metric uncertainties of 5.2% and 2.6% in the FUV and 
NUV, respectively (Morrissey et al. 2007), 2% - 5% in 
the optical/near-infrared bands, and 3% in the 2MASS 
and WISE near- and mid-infrared photometric bands. 

4.2. A N on- Parametric Estimate of the Stellar Mass 
Function 

We build the SMF using the non-parametric 1/U ma x 
estimator widely used in analyses of the galaxy lumi- 
nosity function (see the recent review by Johnston 2011, 
and references therein) . We use this technique in favor of 
complementary parametric methods because accurately 
fitting the observed SMF at z = — 1 requires addi- 
tional free parameters at the low-mass end (Baldry et al. 
2008; Drory et al. 2009; Peng et al. 2010), which cannot 
be reliably constrained at higher redshift with PRIMUS. 
Nevertheless, when calculating number and stellar mass 
densities, we do rely on either a single standard Schcchtcr 
(1976) function, or a double Schechter function (see, e.g., 
Baldry et al. 2008), depending on which model is a bet- 

35 http:/ /www. ucolick.org/~cconroy/FSPS. html 



ter fit to the data, to extrapolate the observed binned 
SMF over small intervals of stellar mass as needed. We 
will also occasionally refer to the 'knee' of the SMF as 
Ai*, which marks the stellar mass above which the SMF 
declines exponentially. 
The differential, non-parametric SMF is given by 

N 

$(logX)A(logX) = ^-^, (3) 

i— 1 

where V max is the maximum cosmological volume within 
which each galaxy i could have been observed given 
the apparent magnitude limits of the survey, uii is 
the statistical weight for each object (see Section 3.1), 
$(log M.) x A(log M.) is the number of galaxies (./V) per 
unit volume with stellar masses in the range logA'J — > 
log M + A(log M), and M has units of M & . 

To estimate V max for each galaxy we first use 
K-correct 36 (v4.2; Blanton & Roweis 2007) to derive 
the redshift-dependent JsT-correction, K(z), from the ob- 
served SED. We use K-correct because of its speed and 
convenience, although we obtain similar results if we in- 
stead use the best-fitting iSEDf it model (Section 4.1). 
Next, we write the apparent magnitude, m, of a galaxy 
of absolute magnitude M as 

m = M + DM + K- Qz, (4) 

where DM(z) is the distance modulus (Hogg 1999) and 
Q is a (constant) luminosity evolution term which we 
discuss below. Given the observed apparent magnitude, 
m Q bj, of an object at redshift z bj, we can use equa- 
tion (4) to write Am(z), the change in the apparent 
magnitude of the object function of redshift, as 

Am(z) = m obj -m(z) = DM(z obj ) - DM(z) 

+ K(z ohj ) - K{z) - Q(z ohj - z). (5) 

By definition, z max (z m \ n ) is the redshift at which A?n 
equals the faint (bright) apparent magnitude limit of 
the survey. Once z max and z m i n have been computed 
for each galaxy with a measured redshift in the interval 
flower, z U pper], Kiax can be computed using 

where z\ = max (z min , zi ower ), z 2 = min (z max , z U p pcr ), Vt 
is the solid angle of the survey, and V c is the comoving 
volume (Hogg 1999). 

The coefficient Q in equations (4) and (5) allows us 
to include a simple luminosity evolution model into 
our Umax estimates. Many recent measurements of 
the optical luminosity function have shown that galax- 
ies brighten toward higher redshift (e.g., Blanton et al. 
2003c; Faber et al. 2007; Loveday et al. 2012; Cool et al. 
2012); consequently, higher-redshift galaxies will be ob- 
servable over a larger cosmological volume (i.e., U max will 
be larger) relative to the no-evolution case. Using Q > 
in equation (5) allows us to mimic the observed luminos- 
ity evolution, at least in a statistical sense. 

Of course, our luminosity evolution model is intention- 
ally simplistic and should not be over-interpreted. For 

36 http:/ /howdy.physics. nyu.edu/index.php/Kcorrect 
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Fig. 2. — Stellar mass vs. rcdshift for all galaxies in the five individual PRIMUS fields. The contours enclose 30%, 60%, and 90% of 
the sample, while the grey points show galaxies lying outside the 90% quantile. The solid black, dot-dashed red, and dashed blue lines 
indicate the stellar mass completeness limit in each field, and the black squares, red diamonds, and blue points show the stellar mass 
limit at the center of each of our six adopted redshift bins (see Table 2). In the CDFS and ELAIS-S1 fields the stellar mass limits among 
all, quiescent, and star-forming galaxies are comparable at all rcdshifts because these two samples are limited in both the optical and at 
3.6 fim (see Section 3.1). Meanwhile, in the COSMOS, XMM-SXDS, and XMM-CFHTLS fields we find that the stellar mass limit for 
quiescent (star-forming) galaxies is higher (lower) at all rcdshifts, as expected given their typically larger (lower) stellar mass-to-light ratios. 
Moreover, the stellar mass limit for the all sample tracks low-mass star-forming galaxies at low redshift, and massive quiescent galaxies at 
high redshift. 



example, Blanton et al. (2003c) find that among z sa 0.1 
galaxies Q varies by roughly a factor of two at rest- 
frame optical wavelengths redward of the 4000-A break, 
whereas our simple model assumes that Q is independent 
of wavelength. Nevertheless, the distribution of V c /V max 
for our SDSS- GALEX sample clearly shows the need to 
account for luminosity evolution when estimating Vmax- 
For example, without luminosity evolution (i.e., Q = 0), 
V c /V max for our SDSS- GALEX sample correlates with 
redshift, stellar mass, and other intrinsic galaxy prop- 
erties; adopting Q = 1 — 2 mag z -1 , on the other 
hand, removes these first-order dependencies and results 
in (V c /V max ) ps 0.5, as expected for a homogenous, sta- 
tistically complete sample (Johnston 2011). The distri- 
bution of V c /V max for PRIMUS, on the other hand, is 
largely insensitive to Q because of the fairly narrow red- 
shift bins we adopt (see, e.g., Figure 1). Using Q > in 
place of Q = changes for just ps 10% of PRIMUS 
galaxies in each redshift interval, and has no significant 
effect on our derived SMFs. 

With the preceding discussion in mind, we calculate 
Knax for each galaxy using equation (6) and the solid 
angle and optical apparent magnitude limits listed in Ta- 
ble 1. We adopt Q = 1.6 mag z _1 (Blanton et al. 2003c; 
Cool et al. 2012) when calculating V max based on our op- 
tical magnitude limits; in the CDFS and ELAIS-S1 fields 
we separately estimate V max based on our 3.6 fim flux 



limits assuming Q = 1.2 mag z _1 (Dai et al. 2009), and 
adopt the smaller of the two (optical vs. mid-infrared) 
Knax values. 

Neglecting stellar mass uncertainties (see Appendix B), 
the two principal sources of uncertainty in the SMF arc 
due to sample size (i.e., Poisson uncertainty) and sam- 
ple variance. In the limit N 3> 1, the formal Poisson 
uncertainty is given by 



<T([> 



A{\ogM)\ 



N 

V 2 

i—l max,: 



(7) 



Equation (7) becomes increasingly inaccurate, however, 
as the number of galaxies approaches zero; therefore, we 
calculate the effective number of galaxies in each mass 
bin following Zhu et al. (2009), and use the analytic for- 
mulae of Gehrels (1986) to compute the upper and lower 
statistical uncertainty of the SMF. 

To estimate the uncertainty in the SMF due to sam- 
ple variance we use a standard jackknifc technique. We 
construct the SMF excluding one field at a time, and cal- 
culate the uncertainty in the mean number of galaxies in 
each stellar mass bin due to the field-to-field variations. 
Formally, we estimate the uncertainty, <j 3 cw , in the j th 
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stellar mass bin due to sample variance as 



°^ = \ ^-E^-^')) 2 ' ( 8 ) 

\ k=l 

where the sum extends over all M individual fields, and 
(3> J } is the mean number density of galaxies in that stel- 
lar mass bin measured from all the available data. Note 
that when computing the cumulative number and stellar 
mass densities (see Section 5.2), we integrate each jack- 
knife realization of the SMF, and compute the variance 
in these quantities using the same formalism. 

This jackknifc technique likely underestimates the level 
of sample variance in PRIMUS, because two of the five 
fields (XMM-SXDS and XMM-CFHTLS) are spatially 
adjacent, so there is covariance between the fields due 
to structure on scales larger than the combined field. 
For simplicity, however, we ignore this covariance in our 
analysis, and simply exclude each of the M = 5 fields 
sequentially as if they were independent. Another po- 
tential issue is that the individual fields are not the same 
size (angular area); therefore, one might expect that the 
standard jackknife prefactor yj {M — 1)/M is not strictly 
correct. However, we verified using a Monte Carlo cal- 
culation that the correct prefactor for our sample differs 
by < 5% from the nominal value despite the more than 
factor of two variation in solid angle among our five fields 
(see Table 1). 

Finally, to estimate the level of sample variance in our 
SDSS- G ALEX sample, we divide the sky into a 12 x 9 
rectangular grid, and retain the 27 30 x 60 deg 2 re- 
gions containing at least 1000 galaxies. We then com- 
pute the variance in the SDSS- GALEX SMF using the 
same methodology described above, once again ignoring 
the potential covariance between adjacent subfields. 

4.3. Stellar Mass Completeness Limits 

Before computing the SMF we need to determine the 
stellar mass above which our sample is complete. In a 
magnitude-limited survey such as PRIMUS, the stellar 
mass completeness limit is a function of redshift, the 
apparent magnitude limit of the survey, and the typi- 
cal stellar mass-to-light ratio of galaxies near this flux 
limit. For example, a quiescent galaxy of a given stellar 
mass will preferentially fall below the survey flux limit 
compared to a star-forming galaxy with the same stellar 
mass, because the stellar mass-to- light ratios of quiescent 
galaxies are typically higher. 

We empirically determine the stellar mass complete- 
ness limits of our sample following Pozzetti et al. (2010). 
First, we compute Mum, the stellar mass each galaxy 
would have if its apparent magnitude was equal to the 
survey magnitude limit, log Mn m = log M + 0.4 (to — 
mu m ), where M. is the stellar mass of the galaxy in units 
of Mq, to is the observed apparent magnitude in the 
selection band, and mu m is the corresponding magni- 
tude limit (see Table 1). Next, we construct the cumu- 
lative distribution of Mu m for the 15% faintest galax- 
ies in Az = 0.04 wide bins of redshift, and calculate 
the minimum stellar mass that includes 95% of the ob- 
jects. We use the subset of galaxies near the flux limit 
to account for the fact that in a flux-limited sample the 
lowest-luminosity galaxies tend to have the lowest stellar 



mass-to-light ratios; however, we did verify that using the 
whole sample in each redshift slice changes the derived 
mass limits by < 0.1 dex. Finally, we fit the limiting stel- 
lar mass versus redshift with a quadratic polynomial sep- 
arately for all, star-forming, and quiescent galaxies. We 
evaluate this fit at the center of each of our six adopted 
redshift intervals (see Figure 1), and list the results in 
Table 2. 

In the CDFS and ELAIS-S1 fields we carry out the 
same procedure described above except we consider the 
flux limits in both the i' and R selection bands, respec- 
tively, and in the IRAC 3.6 fim band (see Section 3.1 and 
Table 1). At each redshift we then take the greater of 
the two stellar mass limits implied by the two apparent 
magnitude limits. In detail, in these fields our sample is 
limited by the 3.6 fim flux limit at low redshift, and by 
the optical flux limit at higher redshift, with the transi- 
tion redshift occurring around z ~ 0.6. 

Figure 2 plots stellar mass versus redshift for the galax- 
ies in all five PRIMUS fields. The black squares, red di- 
amonds, and blue points indicate the stellar mass com- 
pleteness limits at the center of each redshift bin for all, 
quiescent, and star-forming galaxies, respectively, and 
the solid black, dot-dashed red, and dashed blue lines 
show the corresponding polynomial fits. As expected, the 
completeness limits for the quiescent galaxies lie above 
the star-forming galaxy limits at all rcdshifts, except in 
our CDFS and ELAIS-S1 fields, which arc flux-limited 
in both the optical and at 3.6 fim. Moreover, the com- 
pleteness limits for all galaxies overlap the star-forming 
galaxy limits at low redshift, and the quiescent galaxy 
limits at high redshift. This shift occurs because at low 
redshift the combined sample is dominated by low-mass 
star-forming galaxies, while at high redshift massive, qui- 
escent galaxies dominate, as we demonstrate in the next 
section. 

Finally, for our SDSS- G ALEX sample we adopt a 
uniform stellar mass limit of 10 9 M@, which is safely 
above the surface brightness and stellar mass-to-light 
ratio completeness limits of the survey (Blanton et al. 
2005a; Baldry et al. 2008). 

5. EVOLUTION OF THE STELLAR MASS 
FUNCTION FROM z = - 1 

We now have all the ingredients needed to compute 
the evolution of the SMF from z = — 1. We begin 
in Section 5.1 by presenting the SMF at z w 0.1 using 
our SDSS- GALEX sample. Next, we combine the SDSS- 
GALEX and PRIMUS samples in Section 5.2 and show 
how the SMF of all, quiescent, and star-forming galaxies 
has evolved since z — 1. The results wc present in this 
section are all based on our fiducial stellar mass estimates 
(see Section 4.1), although in Appendix B we present a 
detailed discussion of how systematic uncertainties in our 
stellar mass estimates affect our conclusions. 

5.1. Stellar Mass Function at z « 0.1 

We begin by presenting in the upper panel of Fig- 
ure 3 the SMF at z « 0.1 based on the entire SDSS- 
GALEX sample (black squares), and separately for the 
star-forming (blue points) and quiescent (red diamonds) 
galaxy subsamples. Each point represents the comoving 
number density of galaxies in 0.1 dcx wide bins of stellar 
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Fig. 3.— (Upper panel) SDSS-GALEX SMF for all (black 
squares), quiescent (red diamonds), and star- forming (blue points) 
galaxies at z rj 0.1. (Lower panel) Fraction of quiescent galax- 
ies as a function of stellar mass. The massive end of the SMF 
is overwhelmingly comprised of quiescent galaxies, while below 
M ~ 3 X 10 Mq star-forming galaxies increasingly dominate 
the global galaxy population. Quantitatively, the fraction of quies- 
cent galaxies ranges from ~ 25% around ~3x 10 9 A4q to ~ 95% 
around ~3x 10 11 Mq. 

mass, and the vertical error bars indicate the quadrature 
sum of the Poisson and sample variance uncertainties in 
each stellar mass bin. We tabulate the SMF for each 
sample in Table 3. The lower panel of this figure shows 
the variation in the fraction of quiescent galaxies with 
stellar mass. 

Figure 3 conveys several striking (albeit well-known) 
results. First, the massive end of the SMF is al- 
most entirely comprised of quiescent galaxies, while star- 
forming galaxies vastly outnumber quiescent galaxies at 
the low-mass end (see, e.g., Blanton & Moustakas 2009, 
and references therein). Above ~ 2 x 10 11 Mq more 
than ~ 90% of galaxies are quiescent, whereas below 
~ 10 10 Mq star- forming galaxies outnumber quiescent 
galaxies by more than a factor of three. The stellar mass 
at which each population begins to outnumber the other 
is M ~ 3 x 10 10 Mq, in good agreement with previous 
studies (e.g., Bell et al. 2003; Kauffmann et al. 2003b; 
Baldry et al. 2004). Integrating the observed distribu- 
tions above M = 10 9 Mq yields a total stellar mass 
density of 2.36 x 10 8 Mq Mpc -3 , of which approxi- 
mately 60% resides in quiescent galaxies. 37 For com- 
parison, Baldry et al. (2004) find that 54% - 60% of the 
stellar mass density at z ~ 0.1 is in red, early- type galax- 
ies, where the precise result depends on the method used 
to derive stellar masses. 

We compare our results to previously published mea- 

37 Note that galaxies with M < 10 9 A4q contribute a negligible 
amount to the overall stellar mass budget of the nearby Universe 
(see also Brinchmann ct al. 2004). 



suremcnts of the local SMF in Figure 4, adjusting 
where necessary for differences in the adopted IMF and 
cosmological parameters. We plot the SDSS-GALEX 
SMF using filled black squares, and the results from 
Cole et al. (2001), Bell et al. (2003), Li & White (2009), 
and Baldry et al. (2012) using orange diamonds, red 
circles, green triangles, and blue squares, respectively. 
Overall, our results agree reasonably well with these 
studies, although there are some notable differences. The 
agreement between our SMF and the recent measure- 
ment by Baldry et al. (2012), who analyzed a sample of 
~ 10 5 galaxies at z < 0.06 over 143 deg 2 with spectro- 
scopic redshifts from the SDSS and GAMA (Driver ct al. 
2011) surveys, is especially good. Unfortunately, the 
Baldry et al. (2012) sample included too few galaxies 
with stellar masses M >3x 10 11 Mq for them to reli- 
ably measure the massive end of the SMF. 

Compared to Li & White (2009), the exponential tail 
of our SMF falls off less steeply, which is somewhat sur- 
prising given that they analyzed a comparably large sam- 
ple of SDSS galaxies. However, Bernardi et al. (2010) 
argue that Li & White likely underestimated the stel- 
lar masses of the most massive galaxies in their sample 
for two reasons: first, Li & White used Petrosian mag- 
nitudes, which are known to underestimate the fluxes of 
galaxies with extended surface brightness profiles such as 
the spheroidal galaxies that dominate the massive end of 
the SMF (see also Blanton et al. 2011); and second, Li 
& White derived stellar masses using the standard set 
of K-correct basis templates (Blanton & Roweis 2007), 
which can underestimate the stellar masses of massive 
early-type galaxies dominated by very old stellar popu- 
lations (see Bernardi et al. 2010 and Appendix B). 

Finally, Figure 4 shows that the SMF measured by 
Bell et al. (2003) lies systematically above our SMF at 
all stellar masses. Bell et al. constructed their SMF from 
a sample of ~ 7000 galaxies distributed over ~ 400 deg 2 
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for all galaxies against previous determinations from the litera- 
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with some notable differences (see Section 5.1). 
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Fig. 5.— Comparison of the SDSS-GALEX SMF of star- 
forming (top) and quiescent (bottom) galaxies against the SMFs 
of optically-selected blue-cloud and red-sequence galaxies from 
Bell et al. (2003) and Baldry et al. (2012), respectively. We ver- 
ified that when we divide the galaxy population using the u — r 
vs. M r color-magnitude diagram we obtain excellent agreement 
with the corresponding SMFs from Baldry et al.. Therefore, we 
attribute the apparent differences between the SMFs to contami- 
nation of the red sequence by dusty star-forming galaxies. 



previous analyses have used the optical color-magnitude 
diagram to identify "quiescent" and "star-forming" 
galaxies according to whether they lie on the red sequence 
or the blue cloud (e.g., Bell et al. 2003; Taylor et al. 2009; 
Baldry et al. 2012); however, the optical red sequence is 
known to contain both bona fide quiescent galaxies and 
dust-obscured star-forming galaxies (e.g., Brand et al. 
2009; Mailer et al. 2009; Zhu et al. 2011). By contrast, 
we identify quiescent and star-forming galaxies accord- 
ing to whether they lie on or below the star-forming 
sequence, which allows us to select a purer sample of 
quiescent galaxies (Salim et al. 2007; Schiminovich et al. 
2007). Therefore, we anticipate that the overall normal- 
ization of our quiescent-galaxy SMF will be lower relative 
to these previous studies. Moreover, the differences are 
likely to be stellar mass-dependent because the level of 
star formation activity and amount of dust attenuation 
vary systematically with stellar mass (Brinchmann et al. 
2004; Elbaz et al. 2007; Garn & Best 2010). 

With the preceding ideas in mind, in Figure 5 we com- 
pare our SMFs for star-forming and quiescent galaxies 
against the SMFs published by Bell et al. (2003) and 
Baldry et al. (2012). Bell et al. used the g — r color- 
magnitude diagram to identify red-sequence and blue- 
cloud galaxies, while Baldry et al. leveraged the bi- 
modality in u — r color. As expected, our SMF of quies- 
cent galaxies agrees reasonably well with the SMFs from 
Bell et al. and Baldry et al. at the massive end, where the 
amount of contamination from dusty starburst galaxies 
is minimal, but is displaced systematically below their 
SMFs below ~ 10 11 M.q- As anticipated above, the 
reason for these differences is likely because the g — r 
and u — r red sequences at intermediate mass contain 
an admixture of both quiescent and dust-obscured star- 
forming galaxies. In our analysis, these galaxies are (cor- 
rectly) assigned to the star-forming galaxy SMF, as illus- 
trated in the upper panel of Figure 5. We verified this 
interpretation by constructing the SMFs for red-sequence 
and blue-cloud galaxies in our SDSS-GALEX sample us- 
ing the same u — r versus M r optical color-magnitude 
diagram as Baldry et al., and found outstanding agree- 
ment. 



in the SDSS Early Data Release (EDR) survey area 
(Stoughton et al. 2002). However, this area of the sky is 
now known to contain one of the largest structures ever 
mapped, the SDSS Great Wall at z = 0.078 (Gott et al. 
2005). Therefore, one possibility for the origin of the 
discrepancy is that the Bell et al. SMF may be more 
affected by large-scale structure (i.e., sample variance) 
than originally estimated. Alternatively, a strong color- 
dependent difference in stellar mass-to-light ratio that is 
negligible for massive quiescent galaxies and ~ 0.3 dex 
(factor of ~ 2) for low-mass star-forming galaxies could 
also explain the observed discrepancy (E. F. Bell 2012, 
private communication). Whatever the reason, a prac- 
tical consequence of this result is that previous studies 
which relied on the Bell et al. SMF, and to a lesser 
degree the Cole et al. SMF as their low-redshift an- 
chor, may have overestimated the amount of number and 
stellar mass density evolution (e.g., Brammer et al. 2011; 
Mortlock et al. 2011). 

Next, we compare our SMFs of quiescent and star- 
forming galaxies against previous measurements. Most 



5.2. Evolution of the Stellar Mass Function 
Since z w 1 

In the previous section we established the SMF at 
0.1 using our SDSS-GALEX sample. Here, we measure 
how the stellar mass distribution of all, quiescent, and 
star- forming galaxies has changed since 

5.2.1. All Galaxies 

We begin by comparing the SMFs in the five individ- 
ual PRIMUS fields. In Figure 6 we plot the SMFs in 
0.15 dex wide bins of stellar mass divided into six red- 
shift bins from z = 0.2 — 1.0 centered on (z) = 0.25, 0.35, 
0.45, 0.575, 0.725, and 0.9. We choose these redshift bins 
because they correspond to roughly equal ~ 0.9 Gyr in- 
tervals of cosmic time. For clarity we only plot the por- 
tion of each SMF above our stellar mass completeness 
limit in each field (see Section 4.3). 

We find good overall agreement among the individ- 
ual SMFs, modulo expected deviations due to sample 
variance. In Figure 7 we illustrate the effects of large- 
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Fig. 6. — Comparison of the SMFs in each of the five individual PRIMUS fields in six redshift bins from z = 0.2 — f .0. The error bars 
reflect the statistical (Poisson) uncertainty in each stellar mass bin, and for clarity we only show the SMF in each field above our stellar 
mass completeness limit (see Section 4.3). This comparison demonstrates the overall consistency of the SMFs across the five fields, modulo 
expected differences due to sample variance. 
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Fig. 7. — Projected surface density of galaxies vs. redshift in 
the five individual PRIMUS fields as indicated in the legend, and 
for the combined sample (thick black histogram). The vertical 
dotted lines indicate the boundaries of the six redshift bins we have 
adopted. The variation in the number density of galaxies in each 
field and redshift interval bin due to large-scale structure (sample 
variance) is striking. The largest overdensities are in the COSMOS 
field atzs: 0.35, 0.7, and ss 0.85; in the XMM-CFHTLS field 
at z « 0.45, Ri 0.6, and RJ 0.75; and in our CDFS field at z M 0.7. 
Note that the overall surface density in our CDFS and ELAIS-S1 
fields is lower at all redshifts because of the additional 3.6 (im flux 
limit imposed in these fields (see Section 3.1). 



scale structure explicitly by plotting the differential sur- 
face density of galaxies brighter than i ps 23 versus red- 
shift in each of our five fields, and for our combined 
PRIMUS sample. 38 We find significant overdensities in 
the COSMOS field at z 0.35, w 0.7, and w 0.85 (see 
also Lilly et al. 2009; Kovac et al. 2010); in our XMM- 
CFHTLS field at z « 0.45, « 0.6, and w 0.75; and in 
our CDFS field at z « 0.7. By constructing the area- 
weighted average of all five fields (thick black histogram) 
we are able to reduce these field-to-ficld variations sig- 
nificantly, although the effects of sample variance are 
still apparent. Indeed, in this and subsequent sections 
we show that even with five independent fields cover- 
ing w 5.5 deg 2 , sample variance frequently limits the 
precision with which we can constrain the evolution of 
the SMF since z = 1. This conclusion is particularly 
sobering when one considers that all previous analyses 
of the SMF at intermediate redshift which utilized spec- 
troscopic redshifts have been based on samples covering 
at most 1 — 2 deg 2 . In any case, in the remainder of this 
paper we analyze the SMF constructed from the area- 
weighted average of all five fields, and use the jackknife 
technique described in Section 4.2 to empirically estimate 
the uncertainty in the SMF due to sample variance. 

In Figure 8 we plot the SMF of all galaxies in seven 
redshift bins from 2 = — 1 using our combined SDSS- 
GALEX and PRIMUS samples (see Table 4). The black 
squares reflect the comoving number density of galaxies 

38 Recall that an [3.6] < 21 flux cut was also applied to the 
CDFS and ELAIS-S1 samples (see Section 3.1). 




Fig. 8. — Evolution of the SMF from 2 = — 1. The black squares show the comoving number density of galaxies in 0.1 dex wide bins 
of stellar mass based on our SDSS-GALEX (upper-left panel) and PRIMUS samples (subsequent six panels), respectively. Filled (open) 
squares indicate stellar mass bins above (below) the stellar mass completeness limit at the center of each redshift bin. The shaded tan region 
in each panel reflects the quadrature sum of the Poisson and sample variance uncertainties in the SMF, and the solid curve, reproduced 
in every panel for reference, shows the SDSS-GALEX SMF. We find that the SMF for the ensemble population of galaxies has evolved 
remarkably little over the range of redshifts and stellar masses where PRIMUS is complete. 



in 0.1 dcx wide bins of stellar mass, and the filled (open) 
symbols correspond to stellar mass bins above (below) 
our completeness limit in each redshift interval. The tan 
shaded region shows the total (Poisson plus sample vari- 
ance) uncertainty in the SMF, and the solid curve in 
each panel shows the SDSS-GALEX SMF as the null- 
evolution hypothesis. 

Examining Figure 8, we find strikingly little evolution 
in the SMF for the global galaxy population since zfsl, 
at least over the range of stellar masses where PRIMUS 
is complete. In every redshift bin the observed SMF lies 
very close to the local SMF, a result which we quantify 
below. Given the expected stellar mass growth of galax- 
ies due to star formation (see, e.g., Figure 1) and galaxy 
mergers, this result may at first appear surprising. How- 
ever, in Section 5.2.2, we show that the lack of significant 
evolution in the SMF for the ensemble galaxy population 
is a consequence of how the SMFs of the star- forming and 
quiescent galaxy populations separately evolve. More- 
over, in Section 6.1 we show that the relative lack of evo- 
lution in the global SMF, especially at the massive end, 
suggests that mergers do not appear to play a significant 
role for the stellar mass growth of galaxies at z < 1. 

By integrating the SMF above various stellar mass 
thresholds we can quantify the observed (lack of) evo- 
lution in the SMF, and look for evidence of mass assem- 
bly downsizing within the global galaxy population (see 
Section 1). In Figure 9 we plot versus redshift the cu- 



mulative number density of galaxies with stellar masses 
greater than 10 9 5 , 10 10 , 10 10 5 , and 10 11 M Q . We focus 
here on the number density evolution, although the evo- 
lution in stellar mass density leads to the same basic con- 
clusions (see Table 5). We integrate the observed SMF 
directly, but exclude stellar mass bins containing fewer 
than three galaxies where the SMF is noisiest. We use the 
best-fitting Schechter model to extrapolate the observed 
SMF as needed over small intervals of stellar mass. We 
emphasize, however, that these mo del- dependent correc- 
tions typically modify the measured number densities by 
< 0.02 dex, and therefore potential errors in the extrap- 
olations do not affect any of our conclusions. The solid 
black squares in Figure 9 show the mean number density, 
while the vertical error bars indicate the Poisson uncer- 
tainty; the thin black boxes around each point indicate 
the quadrature sum of the Poisson and sample variance 
uncertainties in the vertical direction, and the redshift 
bin width in the horizontal direction. This graphical 
representation shows that sample variance uncertainties 
are frequently comparable to or larger than the Pois- 
son uncertainties. Finally, symbols with upward-pointing 
arrows represent lower limits, and the grey shaded re- 
gion shows for reference the comoving number density of 
galaxies at z « 0.1, to illustrate the case of no evolution. 

Figure 9 shows that the cumulative number density of 
galaxies above all four stellar mass thresholds does not 
appear to change significantly over the range of redshifts 
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Fig. 9. — Evolution of the cumulative comoving space density of galaxies more massive than (top-left) 10 Mq, (top-right) 10 Mq, 
(bottom-left) 10 10 ' 5 Mq, and (bottom-right) 10 11 Mq from 2 = — 1 based on the SMFs presented in Figure 8. The error bars reflect 
the Poisson uncertainty on each number density measurement, and the thin black boxes reflect the quadrature sum of the Poisson and 
sample variance uncertainties in the vertical direction, and the redshift bin width in the horizontal direction. We designate lower limits 
using upward-pointing arrows. The grey shaded region in each panel shows — as the null-evolution hypothesis — the number density of 
galaxies at z 0.1 based on our SDSS- GALEX sample. We find that the cumulative number density of galaxies with M > 10 9 ' 5 Mq and 
M > 10 10 Mq has increased by 15% ± 10% and 21% ± 19% since z = 0.4 and z = 0.6, respectively, while the cumulative space density of 
M > 10 10 ' 5 Mq and M > 10 11 Mq galaxies has changed by just 4.7% ± 12% and 11% ± 17% since z = 0.8 and z = 1, respectively. 



where PRIMUS is complete. To quantify this result, 
we fit a power-law function of redshift, n cx (1 + z) 7 , 
to the measured densities, excluding lower limits. We 
scale the formal statistical uncertainties by the 
square root of the x 2 statistic divided by the number 
of degrees-of-frecdom, in order to be able to intercom- 
pare the significance of the evolutionary trends across all 
four stellar mass thresholds. We find 7 = —0.43 ± 0.3, 
-0.40 ± 0.4, 0.17 ± 0.4, and 0.32 ± 0.4 for the evolution 
in the cumulative number of galaxies more massive than 
10 9 - 5 , 10 10 , 10 10 - 5 , and 10 11 Mq, respectively. Above 
the two highest stellar mass thresholds, the measured 
number densities at z w 0.7, and to a lesser extent at 
z w 0.9, are clearly affected by the large-scale overden- 
sities in several of the PRIMUS fields (see Figure 7). 
Therefore, excluding the z ~ 0.7 redshift bin, we obtain 
7 = -0.08±0.2 and 0.16±0.3 above 10 la5 and 10 11 Mq, 
respectively. Stated another way, the cumulative number 
of M > 10 9 - 5 Mq and M > 10 10 Mq galaxies has in- 
creased by 15% ± 10% and 21% ± 19% since z = 0.4 and 
z = 0.6, respectively. Meanwhile, the cumulative space 



density of M > 10 10 5 Mq and M > 10 11 Mq galax- 
ies has remained relatively constant, changing by just 
4.7% ± 12% and 11% ± 17% since z = 0.8 and z = 1, 
respectively. 

Thus, while we find hints of mass assembly 
downsizing — a more rapid increase in the number of 
lower-mass galaxies toward low redshift — the trends are 
only marginally significant. By contrast, previous studies 
have found much stronger evidence for downsizing within 
the global galaxy population (see, e.g., Fontana et al. 
2006; Perez-Gonzalez et al. 2008; Pozzetti et al. 2007, 
2010). For example, Pozzetti et al. (2010) report a 
32% ± 6% increase in the cumulative number of M > 
10 9 ' 5 Mq galaxies since z = 0.44, and no statistically 
significant change (7% ± 17%) in the number density of 
galaxies more massive than 10 11 Mq since 2=1 based 
on an analysis of the COSMOS field. While our findings 
are qualitatively consistent with these studies — many of 
which pushed further down the stellar mass function at 
higher redshift, and therefore had a larger lever arm with 
which to detect downsizing — we also find that sample 




Fig. 10. — Evolution of the SMFs of quiescent (dark red diamonds and light red hatched region) and star-forming (dark blue points 
and light blue shaded region) galaxies from z = — 1. Filled (open) symbols correspond to stellar mass bins above (below) our stellar 
mass completeness limit in each redshift interval. The dashed blue and solid red curves are the SDSS-GALEX star-forming and quiescent- 
galaxy SMFs (upper-left panel), and have been reproduced in every panel for reference. We find a significant increase in the number of 
intermediate-mass (~ 10 10 Mq) quiescent galaxies toward lower redshift, but essentially no change in the SMF of quiescent galaxies above 
~ 10 11 Mq. Meanwhile, the SMF of star-forming galaxies is largely invariant below ~ 10 11 Mq (at least where PRIMUS is complete), 
but exhibits significant evolution above ~ 10 Mq by shifting toward lower stellar mass at fixed number density with decreasing redshift. 



variance can wash out the significance of the observed 
trends (see also Fontanot et al. 2009). In any case, we 
will show in Section 5.2.2 that by measuring the change 
in the number density of galaxies within fixed bins of stel- 
lar mass (as opposed to using stellar mass thresholds) the 
signatures of downsizing will become more apparent. 

5.2.2. Quiescent and Star-Forming Galaxies 

In the previous section we measured the evolution of 
the SMF for the ensemble population of galaxies from 
z = — 1. We found only a marginally significant in- 
crease (16% ± 9%) in the cumulative number density 
of all M > 10 9 - 5 Mq galaxies since z 0.6, and 
very little change (8% ± 10%) in the space-density of 
M > 10 10 ' 5 .M© galaxies since 2 » 1. Here, we investi- 
gate these results in more detail by dividing our sample 
into quiescent and star- forming galaxies based on the cri- 
teria defined in Section 3.2. 

In Figure 10 we plot the SMFs of quiescent and 
star-forming galaxies in seven redshift bins from z = 
— 1. Filled (open) symbols indicate stellar mass bins 
above (below) our completeness limit at the center of 
each redshift bin, and the solid red and dashed blue 
curves, reproduced in every panel for reference, show 
the SDSS- GALEX quiescent and star-forming galaxy 
SMFs, respectively. In Section 5.3 we compare our type- 



dependent SMFs against published measurements and 
show that they are broadly consistent with previous stud- 
ies. 

Figure 10 shows that the strong stellar mass depen- 
dence of galaxy bimodality observed among local galax- 
ies persists over the full range of redshifts probed by 
PRIMUS. In other words, we find that quiescent galax- 
ies dominate the massive end (> 10 11 Mq) of the SMF, 
and star-forming galaxies dominate among intermediate- 
mass (~ 10 10 Mq) galaxies at all redshifts from z = 
— 1. However, we also observe several striking evolu- 
tionary trends. Among quiescent galaxies, the number 
of intermediate-mass galaxies increases dramatically to- 
ward the current epoch, while the massive end of the 
SMF remains remarkably fixed. Meanwhile, the largest 
changes in the SMF of star-forming galaxies occur at 
the massive end. We find a perceptible shift in the 
star-forming galaxy SMF toward lower mass at fixed 
number density with decreasing redshift, while the low- 
mass end of the SMF remains relatively constant over 
the whole range of stellar masses and redshifts where 
our sample is complete. The so-called transition mass — 
the stellar mass at which the quiescent and star-forming 
galaxy SMFs cross — evolves roughly ascx (1 + z) 1 - 5 , from 
- 3 x 10 10 Mq at z w 0.1 to - 7 x 10 10 Mq at 
z sj 0.9, which agrees reasonably well with previous mea- 




Fig. 11. — Evolution of the SMFs of (left) star-forming and (right) quiescent galaxies from z = — 1 based on the data presented in 
Figure 10. We use progressively lighter shades of blue to show the evolution of the star-forming galaxy SMF, and shades of orange to show 
how the SMF of quiescent galaxies has evolved. The black shaded region in each panel shows the corresponding SDSS- GALEX SMF. 



surements (e.g., Bundy et al. 2006; Vergani et al. 2008; 
Pozzetti et al. 2010). It is not clear, however, that the 
transition mass has any physical interpretation, as Fig- 
ure 10 shows that its evolution is entirely driven by the 
rise in the number of intermediate-mass quiescent galax- 
ies (e.g., Borch et al. 2006; Cattaneo et al. 2008). 

Another way to visualize these results is with Fig- 
ure 11, which shows the individual SMFs from all seven 
rcdshift bins on top of one another. For clarity, we only 
plot each SMF above our stellar mass completeness limit, 
and we only show stellar mass bins containing three or 
more galaxies. In the left panel we use progressively 
lighter shades of blue to show the evolution of the star- 
forming galaxy SMF, and in the right panel we show 
the evolution of the quiescent-galaxy SMF using progres- 
sively lighter shades of orange. The black shaded region 
shows the corresponding SDSS- GALEX SMF, which we 
plot on top so that the changes in the SMF with red- 
shift can be more easily evaluated. This figure clearly 
shows the significant steepening of the low-mass end of 
the SMF of quiescent galaxies toward lower redshift, and 
the simultaneous decline in the number of massive star- 
forming galaxies. 

In Figure 12 we quantify the observed evolution by 
plotting the integrated number density of galaxies mea- 
sured in four 0.5 dex wide intervals of stellar mass be- 
tween 10 9 - 5 - 10 u ' s M Q . As in Section 5.2.1, we cal- 
culate the number density by integrating the observed 
SMFs directly, excluding stellar mass bins with fewer 
than three galaxies, and use our Schechter model fits 
to extrapolate over small ranges of stellar mass. We plot 
the evolution in the number density of all, quiescent, and 
star-forming galaxies using black squares, red diamonds, 
and blue points, respectively, and indicate lower limits 



on the number density in redshift bins where our SMF is 
partially incomplete using upward-pointing arrows. The 
error bars reflect the quadrature sum of the Poisson and 
sample variance uncertainties. We list the derived num- 
ber and stellar mass densities in Table 6. 

We quantify the observed trends by fitting a power-law 
function of redshift to the measured number densities, 
given by 

n(z) = n (l + z)\ (9) 

In addition, we model the evolution in the stellar mass 
density, p{z), as 

p(z) = p (l + zf. (10) 

In detail, we fit the data in log (n) — log (1 + z) and 
log (p) — log (1 + z) space using weighted linear least- 
squares minimization, and we only fit over the range of 
rcdshifts where our measurements are complete (i.e., we 
ignore lower limits). We also exclude from the fits our 
measurements at z w 0.7 because of the above- average 
overdensity of galaxies in this redshift bin (see Figure 7) . 
We emphasize, however, that including this redshift bin 
would only strengthen our claim of minimal evolution in 
the number density of massive galaxies; in other words, 
excluding the z ss 0.7 measurements is a conservative 
choice. The solid black, dot-dashed red, and dashed blue 
lines in Figure 12 show the results of fitting the num- 
ber density of all, quiescent, and star-forming galaxies, 
respectively, and the corresponding shaded regions show 
the la range of power-law fits drawn from the full covari- 
ance matrix. Table 7 lists the best-fitting coefficients and 
uncertainties, where the uncertainties have been rescaled 
as in Section 5.2.1 such that xi = 1- 
Figure 12 synthesizes nearly all the key results of this 
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Fig. 12. — Evolution in the number density of all (black squares), quiescent (red diamonds), and star-forming (blue points) galaxies in 
four 0.5 dex wide intervals of stellar mass ranging from 10 9 - 5 - 10 10 Mq in the upper-left panel, to 10 11 — 10 A4q in the lower-right 
panel. The error bar on each measurement is due to the quadrature sum of the Poisson and sample variance uncertainties in each redshift 
interval; we denote lower limits on the number density using upward-pointing arrows. The solid black, dot-dashed red, and dashed blue 
lines show weighted linear least-squares fits to the data, and the corresponding shaded regions show the la range of power-law fits drawn 
from the full covariance matrix. We find a factor of ~ 2 — 3 increase in the number density of ~ 10 9 ' 5 — 10 8 M<3 quiescent galaxies 
since z « 0.5, and a remarkably little change (—8% ± 10%) in the space density of star-forming galaxies over the same range of stellar 
mass and redshift. Between 10 10 ' 5 — 10 11 A4q, the space density of quiescent galaxies increases by 58% ± 9%, while the number density 
of star-forming galaxies declines by —13% ± 23%. Meanwhile, above 10 11 A4q we find a steep 54% ± 7% decline in the number of massive 
star-forming galaxies since z rj 1, and a small increase (22% ± 12%) in the space-density of comparably massive quiescent galaxies. The 
distinct evolutionary trends exhibited by star-forming and quiescent galaxies conspire to keep the number density of all galaxies relatively 
constant over the range of stellar masses and redshifts probed by PRIMUS. 



section, and conveys many of the core conclusions of 
this paper. First, we find that the number density of 
10 95 — 10 10 M.q quiescent galaxies increases significantly 
toward lower redshift, by a factor of 3.2 ± 0.5 since 
z = 0.4, whereas the number density of star-forming 
galaxies decreases marginally, by —10% ± 15% over the 
same redshift range. Meanwhile, the number density of 
10 10 — 10 10 5 A4q quiescent galaxies increases by a factor 
of 2.2±0.4 since z = 0.6, while the number of comparably 
massive star-forming galaxies changes by — 4%±15%. Fi- 
nally, in the 10 10 ' 5 — 10 11 A4q stellar masses bin we find a 
58%±9% increase in the space density of quiescent galax- 
ies since z = 0.8, and a — 13%±23% decrease in the num- 
ber of star-forming galaxies over the same redshift range. 
Thus, we find remarkably little change (—8% ± 10%) in 
the number density of 10 9 5 — 10 10 5 A4q star-forming 
galaxies over the full range of redshifts where PRIMUS 
is complete, and a gradual, but significant buildup in 
the population of quiescent galaxies toward low redshift. 



Moreover, we find that the rate at which the quiescent 
galaxy population builds up toward low redshift increases 
steeply with decreasing stellar mass. 

Among the most massive galaxies in our sample, how- 
ever, Figure 12 reveals a striking inversion of the trends 
seen at lower stellar mass. Between 10 11 — 10 115 M.q 
the number density of quiescent galaxies increases by 
22% ± 12% since z ps 1, while the number of massive 
star-forming galaxies declines by 54%±7% over the same 
redshift range. The reason this destruction of massive 
star-forming galaxies (and presumed transformation into 
quiescent systems) does not significantly affect the space 
density of massive quiescent galaxies is because quiescent 
galaxies vastly outnumber star-forming galaxies above 
~ 10 11 A4q at all redshifts from z = — 1. For example, 
at z = 1 quiescent galaxies outnumber 10 11 — 10 115 Mq 
star- forming galaxies by w 0.13 dex (« 35%); therefore, 
the w 55% decline in the number of massive star-forming 
galaxies can easily be subsumed into the quiescent galaxy 
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Fig. 13. — Comparison of the SMFs for (top row) quiescent and (bottom row) star-forming galaxies in four redshift bins from z = 0.2 — 1.0 
against previously measurements SMFs assembled from the literature. We find that our SMFs are generally consistent with these previous 
studies. 



population by the current epoch. 

Reviewing Figures 10—12, it is now clear why the 
SMF for the global population of galaxies (Figures 8 
and 9) evolves so little since z ps 1 over the range of 
stellar masses where PRIMUS is complete. Between 
lO 9 ' 5 - 10 10 - 5 M Q the SMF is dominated by star- forming 
galaxies, whose number density remains relatively con- 
stant. Meanwhile, among Ai > 10 10,5 A4@ galaxies 
the SMF becomes increasingly dominated by quiescent 
galaxies, whose number density also remains roughly 
constant with redshift. Consequently, the bimodal na- 
ture of the galaxy SMF combined with the distinct evo- 
lutionary trends exhibited by star-forming and quiescent 
galaxies conspire to keep the SMF for the global pop- 
ulation of galaxies from changing significantly at these 
rcdshifts. 

In Section 5.2.1 we found hints of differential evolu- 
tion in the SMF of all galaxies based on stellar mass- 
threshold samples, but the results were not very signifi- 
cant. Do we find more significant evidence for mass as- 
sembly downsizing based on the number densities derived 
within fixed-interval bins of stellar mass? Our power-law 
fits to the black squares in Figure 12 (see also Table 7) 
indicate a 28% ± 11% increase in the space density of all 
10 9 - 5 - 10 10 M Q galaxies since z = 0.4, and a 35% ± 14% 
increase in the number density of 10 10 — 10 10 5 A4q galax- 
ies since z = 0.6. By contrast, among 10 10 5 — 10 11 A4q 
galaxies the number density increases by 4% ± 9% since 



z = 0.8, and declines by -10%±9% for 10 u -10 u ' s M Q 
galaxies since z = 1. Thus, we do find evidence for 
mass assembly downsizing — a continued buildup of the 
low- and intermediate-mass galaxy population toward 
low redshift, and no significant changes in the space den- 
sity of massive galaxies — within the global galaxy pop- 
ulation. However, with the benefit of hindsight we now 
see that these relatively subtle evolutionary trends are 
being driven entirely by the much more significant evo- 
lutionary trends separately exhibited by the population 
of quiescent and star-forming galaxies. 

To summarize, we have shown that the evolution of 
the SMFs of both quiescent and star- forming galaxies de- 
pends sensitively on stellar mass. Above M. ~ 10 11 A4q 
quiescent galaxies dominate the galaxy population at 
all redshifts, and their number density changes rela- 
tively little since z = 1, whereas the number of star- 
forming galaxies declines precipitously toward lower red- 
shift. Between 10 95 — 10 105 M®, on the other hand, 
star-forming galaxies vastly outnumber quiescent galax- 
ies, and their number density changes by just — 8%± 10% 
between z ps 0.6 and z ps 0. Meanwhile, the number 
of 10 95 — 10 10 5 A4q quiescent galaxies increases sig- 
nificantly since z « 0.6 at a rate that accelerates with 
decreasing stellar mass. Finally, the stellar mass range 
between 10 10 5 — 10 11 A4q marks the transition regime 
between the dominance of star-forming galaxies and the 
rise of quiescent galaxies at low mass, and the dominance 
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of quiescent galaxies and the decline of the star-forming 
population at large stellar mass. 

5.3. Comparison With Previous Studies 

A detailed quantitative comparison of our results 
against previous studies is challenging for several rea- 
sons. First, previous studies have used a wide variety 
of techniques to divide the galaxy population into "qui- 
escent" and "star-forming" galaxies; however, many of 
these techniques result in a quiescent population that is 
highly contaminated by dusty star- forming galaxies (e.g., 
Mailer et al. 2009; Zhu et al. 2011), which can severely 
bias the inferred SMF in mass-dependent ways. Sec- 
ond, systematic differences in stellar mass estimates due 
to different prior assumptions and population synthesis 
models can significantly affect the inferred SMF (e.g., 
Marchesini et al. 2009; Kajisawa et al. 2009; see also Ap- 
pendix B) And finally, many previous studies have ne- 
glected the effects of sample variance, and therefore have 
underestimated the statistical uncertainties of their re- 
sults. Nevertheless, we can still perform a rudimen- 
tary comparison of our quiescent and star-forming galaxy 
SMFs against previous measurements assembled from 
the literature. 

To facilitate this comparison, we recompute our SMFs 
using four broader redshift bins with Az = 0.2 (due to 
the typically smaller area and sample size of these anal- 
yses) centered on (z) = 0.3, 0.5, 0.7, and 0.9. In Fig- 
ure 13 we plot the SMFs for (top row) quiescent and (bot- 
tom row) star-forming galaxies from PRIMUS as a tan 
shaded region, reflecting the quadrature sum of the Pois- 
son and sample variance uncertainties in each redshift 
bin. We compare these results to the SMFs published by 
Borch et al. (2006) (red triangles), Drory et al. (2009) 
(purple diamonds), Pozzctti et al. (2010) (blue squares), 
and Ilbert et al. (2010) (orange crosses), accounting for 
differences in the adopted Hubble constant and IMF. For 
reference, the SMFs published by Drory et al. (2009), 
Ilbert et al. (2010), and Pozzetti et al. (2010) are all 
based on the - 2 deg 2 COSMOS field, while Borch et al. 
(2006) analyzed the three COMBO- 17 (Wolf et al. 2003) 
fields, totaling ~ 0.8 deg 2 . 

Examining Figure 13, we find reasonably good agree- 
ment between our SMFs and the literature, modulo ex- 
pected differences due to sample variance in PRIMUS 
(see, e.g., Figure 7) and the other reasons outlined above. 
Among the largest discrepancies are in the (z) = 0.5 
and (z) = 0.7 redshift bins for quiescent galaxies, which 
shows that our SMF agrees with the SMFs derived 
by Borch et al. (2006) and Pozzetti et al. (2010), but 
disagrees noticeably with the Drory et al. (2009) and 
Ilbert et al. (2010) SMFs. We also find a somewhat 
higher number density of ~ 10 10 5 — 10 11 A4q star- 
forming galaxies at (z) — 0.5. Overall, however, we 
conclude that our results are consistent with previous 
measurements of the SMF at intermediate redshift. 

6. DISCUSSION 

We have measured the evolution of the SMF since z m 
1 of quiescent and star- forming galaxies using PRIMUS, 
one of the largest spectroscopic surveys of intcrmcdiatc- 
rcdshift galaxies ever conducted. Our goals have been to 
characterize the stellar mass growth of each population, 



and to measure the rate at which star-forming galaxies 
are being quenched as a function of stellar mass and red- 
shift. Compared to many previous studies, our analysis 
has benefited from a large, statistically complete sample 
of faint galaxies (~ 40, 000 galaxies to i k 23) spread 
across five widely separated fields totaling sa 5.5 deg 2 , 
and a well-defined local SDSS- G ALEX comparison sam- 
ple. With these data, we have been able to study the de- 
tailed evolution of the SMF over a large dynamic range 
of stellar mass and redshift in a consistent way, with well- 
quantified sample variance uncertainties. 

We find that the evolution of the SMFs of both quies- 
cent and star-forming galaxies depends acutely on stel- 
lar mass, but in very different ways (see Figure 12). 
Among quiescent galaxies, the number of intermediate- 
mass (~ 10 10 A4q) galaxies increases by a factor of 
~* 2 — 3 since z « 0.5, but remains approximately con- 
stant for massive (> 10 11 M Q ) galaxies since z m 1. 
By contrast, the most significant evolutionary trends for 
star-forming galaxies occur above ~ 10 11 A4q. Specifi- 
cally, we find no significant change in the number density 
of intermediate-mass galaxies, and a f=a 55% decrease in 
the number of massive star- forming galaxies since zwl. 

These galaxy-type dependent trends conspire rather 
remarkably to make the SMF for the global galaxy pop- 
ulation appear to not have changed significantly since 
2 w 1, at least over the range of stellar masses and red- 
shifts probed by PRIMUS (Figure 8). One implication 
of these results is that an analysis of the global galaxy 
population by itself would yield a highly incomplete view 
of galaxy evolution because it would mask the rich inter- 
play between the coevolution of star-forming and quies- 
cent galaxies. 

In the next two sections we synthesize our results to in- 
vestigate the effect of mergers on the stellar mass growth 
of galaxies, and to quantify the stellar mass dependence 
of star formation quenching from z = — 1. 

6.1. Constraints on the Stellar Mass Growth of 
Galaxies by Mergers 

The stellar mass of an individual galaxy can change 
by forming new stars, or by merging with other galax- 
ies. By accounting for the stellar mass growth of galax- 
ies by in situ star formation, the redshift evolution of 
the SMF (i.e., the derivative of the SMF with respect 
to cosmic time) in principal can be used to constrain 
the growth of galaxies by mergers (e.g., Drory & Alvarez 
2008; Walcher et al. 2008; Conroy & Wechsler 2009; 
Pozzetti et al. 2010). 

Our finding that the SMF of the global galaxy popu- 
lation evolves relatively little between z « 1 and z 
(see Figure 8) suggests that mergers play a subdominant 
role for the stellar mass growth of galaxies at these red- 
shifts. To quantify this result, we use the measured SFR 
of each galaxy (see Section 4.1) to estimate how much 
their stellar mass will increase by in situ star formation. 
Specifically, we compute for each galaxy of a given stellar 
mass, M, a new stellar mass, A4', given by 

M'(z') =M(z) + {l-K) x (tpAt), (11) 

where 1Z is the return fraction, the stellar mass returned 
(assumed instantaneously) to the interstellar medium by 
supcrnovae and stellar winds, ip is the SFR in A4q yr -1 , 
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and At = t{z') — t{z) is the elapsed cosmic time (increase 
in the age of the Universe) in years between redshift z and 
z', where z 1 < z. We adopt 7Z w 0.5, which is appropriate 
for the Chabrier (2003) IMF (Leitner & Kravtsov 2011), 
and make the simplifying assumption that the SFR is 
constant over the time interval At. Although the SFRs of 
most star-forming galaxies at these redshifts are declining 
with decreasing redshift (see, e.g., Figure 1; Noeske et al. 
2007), At in our analysis is sufficiently short (< 1 Gyr) 
for all but the last redshift bin (where At « 1.6 Gyr), 
that incorporating a more detailed star formation his- 
tory into our calculation would not significantly change 
our results. Moreover, recall that our UV-based SFRs 
trace star formation over the last ~ 100 Myr, which is 
reasonably well matched to the At timcscale. 

Using this formalism, we can use the observed SMF 
at redshift z to predict what the SMF will look like at 
a (lower) redshift z' if galaxies grow by star formation 
alone (i.e., assuming mergers do not occur). Note that al- 
though in situ star formation conserves the total number 
of galaxies, the number of galaxies of a given stellar mass 
can increase or decrease because the SFR varies with stel- 
lar mass (see, e.g., Figure 1). After controlling for star 
formation growth, any residual change in the number 
density of galaxies of a given stellar mass must be due 
to merging. We emphasize that our measurement com- 
plements, but is only implicitly related to measurements 
of the major and minor merger rate (see, e.g., Lotz et al. 
2011, and references therein). For example, the tech- 
nique we use only reveals whether mergers have a net ef- 
fect on the SMF, but cannot be used to infer the underly- 
ing stellar mass distribution of galaxies that are merging 
(see Drory & Alvarez 2008 for more details). Another 
caveat regarding this technique is its implicit assump- 
tion that mergers retain all the stellar mass involved in 
the merger, when it is likely that a non-negligible frac- 
tion of that mass is dispersed to large radii to form the 
diffuse stellar component (DSC) of groups and clusters 
(e.g., Murante et al. 2007). In fact, our precise measure- 
ment of the evolution of the SMF above ~ 10 11 Mq 
could be turned around to help constrain this fraction 
in the context of the cosmological growth of dark-matter 
halos (e.g., Monaco et al. 2006; Bchroozi et al. 2012). 

With the preceding discussion in mind, we define the 
fractional merger growth rate, G(M, 2), as 
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(12) 

where n ob8 and n prod are the observed (measured) and 
predicted number density of galaxies at redshift z, re- 
spectively, and At is in Gyr. We divide by the observed 
number density to account for the shape of the SMF; for 
example, a merger-induced absolute increase in number 
density of 10 -3 galaxies Mpc~ 3 Gyr -1 will have a much 
larger fractional effect on the exponential tail of the SMF 
relative to the low-mass end, where such a change would 
be negligible. We derive the number densities in equa- 
tion (12) as in Section 5.2 by numerically integrating 
the observed and predicted SMFs over chosen intervals 
of stellar mass, but use the best-fitting single or double 
Schcchtcr model to extrapolate to lower or higher mass 
as needed. Note that Q can be either positive or negative 
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Fig. 14. — Merger growth rate, G{M,z), the fractional change 
in the number density of galaxies due to mergers after account- 
ing for stellar mass growth by in situ star formation, vs. redshift 
in four 0.5 dcx wide intervals of stellar mass between 10 9 ' 5 Mq 
and 10 11,5 Mq. Note that Q as denned by equation (12) can be 
either negative or positive depending on whether mergers prefer- 
entially destroy or create galaxies of a given stellar mass. We find 
a net change in the number density of 10 9 ' 5 — 10 10,5 Mq galaxies 
of —13% ± 9% Gyr" 1 from z as to z ~ 0.5, and a change of 
-14% ± 18% Gyr" 1 in the number density of 10 10 ' 5 - 10 11 ' 5 Mq 
galaxies since z K, 0.8. In other words, although mergers are almost 
certainly occurring, they do not have a large net effect on the SMF 
over this range of stellar masses and redshifts. 

depending on whether mergers result in a net increase or 
decrease of galaxies of a certain stellar mass. Moreover, 
Q w does not necessarily indicate that mergers are not 
occuring or important; mergers could result in a small 
value of Q if the growth and destruction of galaxies into 
and out of a certain stellar mass bin on average balanced 
one another. 

In Figure 14 we plot Q versus redshift in four intervals 
of stellar mass between 10 9 5 — 10 115 M.q. Qualitatively, 
we find no significant variations of Q with either stellar 
mass or redshift (where PRIMUS is complete). Quanti- 
tatively, we find a net change in the number density of 
10 9.5 _ 1Q io.5 Mq ga i ax ie S f -13% ± 9% Gyr" 1 from 

z w to z w 0.5, and a change of -14% ± 18% Gyr" 1 
in the number density of 10 10 ' 5 — 10 115 M.q galaxies 
since z w 0.8. Although the uncertainties are significant, 
wc conclude, therefore, that mergers do not appear to 
be an important channel for stellar mass growth at late 
cosmic times, even among massive (> 10 11 A4q) galax- 
ies (see also Pozzetti et al. 2010). Although beyond the 
scope of the present study, a detailed comparison of these 
results with theoretical galaxy formation models, which 
generally find that massive galaxies grow much more sub- 
stantially through mergers at z < 1 (e.g., De Lucia et al. 
2006; Cattaneo et al. 2008, 2011), would be very inter- 
esting. 

We extend our analysis further by separately consider- 
ing the stellar mass growth of quiescent and star-forming 
galaxies. One particularly important question we would 
like to address is whether quiescent galaxies grow sig- 
nificantly by dissipationless mergers (also known as dry 
mergers) at z < 1. To investigate this question we mea- 
sure the redshift evolution of Ai c , the stellar mass at 
fixed cumulative number density, n(A4 > Ai c ). To a 
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Fig. 15. — Stellar mass, A4 C , at which the cumulative number 
density, n(M > M c ), equals 10 -3 ' 5 Mpc -3 vs. redshift. The red 
and blue shaded regions correspond to quiescent and star-forming 
galaxies, respectively, and the grey hatched area indicates for ref- 
erence the stellar mass completeness limit for quiescent galaxies; 
the completeness limit for star-forming galaxies extends to much 
lower stellar mass. The constancy of Ai c for quiescent galaxies in- 
dicates negligible growth of this population due to mergers, while 
the decline in JV[ C with decreasing redshift reflects the progressive 
quenching and transformation of massive star-forming galaxies into 
quiescent galaxies from z = — 1. 

good approximation, both star formation and mergers 
will increase M c without changing the number density, 
while the transformation of one galaxy type into an- 
other (e.g., due to star formation quenching) will tend 
to decrease M. c at a given cumulative number density 
(van Dokkum et al. 2010; Brammcr et al. 2011). Be- 
cause the amount of in situ star formation in our sample 
of quiescent galaxies is by construction negligible (see 
Figure 1), an increase in Ai c with decreasing redshift 
can be attributed to dissipationless mergers. 

In Figure 15 wc plot M. c versus redshift for quies- 
cent (hatched red shading) and star-forming (light blue 
shading) galaxies corresponding to n{M > M c ) = 
10~ Mpc~ 3 . We choose this number density thresh- 
old because it samples a significant fraction of the ex- 
ponential tail of the SMF at each redshift, and because 
PRIMUS is complete to both galaxy types all the way to 
z « 1. For reference, the grey shaded region in the lower- 
right corner of Figure 15 shows the stellar mass complete- 
ness limit for quiescent galaxies (the stellar mass limit 
for star-forming galaxies is much lower). In each redshift 
bin, our measurement of M c reflects the quadrature sum 
of both the Poisson and sample variance uncertainties. 

We find no notable change in M. c for quiescent galaxies 
from z = — 1, and a significant decrease for star- forming 
galaxies with decreasing redshift. Modeling the observed 
trends as a power-law function of redshift, M. c oc (l+z) q , 
we find q = —0.09 ± 0.16 for quiescent galaxies, and 
q = 0.60 ±0.17 for star-forming galaxies. In other words, 
M c for star- forming galaxies decreases on average by 
0.18 ± 0.05 dex (51% ± 12%) since 2 = 1, whereas M c 
for quiescent galaxies is constant to within ±0.05 dex 
(±12%) over the same redshift range. 

The decline in M. c for star- forming galaxies is most 
likely due to the progressive transformation of massive 
star-forming galaxies into quiescent, passively evolving 



systems toward low redshift (see Figure 12 and Sec- 
tion 6.2). In essence, M. c decreases toward low redshift 
because we have to integrate further down the SMF to 
count the same total number of galaxies. Meanwhile, 
the constancy of M. c for quiescent galaxies follows di- 
rectly from the lack of evolution at the massive end of the 
SMF. Note that M. c for quiescent galaxies remains ap- 
proximately constant even though star-forming galaxies 
are being quenched because massive star-forming galax- 
ies constitute a very small fraction of the total number 
of massive galaxies at z = — 1 (sec Section 5.2.2). Wc 
conclude, therefore, that most massive, quiescent galax- 
ies are fully assembled by z ~ 1, and do not appear to 
grow significantly by dissipationless mergers over the last 
~ 8 billion years of cosmic time. 

6.2. Buildup of the Quiescent Galaxy Population by 
Star Formation Quenching 

Our analysis of the SMF has revealed significant 
changes in the population of both quiescent and star- 
forming galaxies. The two key results wc focus on in this 
section are the rapid rise in the number of intermediate- 
mass (~ 10 10 Mq) quiescent galaxies since z « 0.5, 
and the steep decline in the population of massive (> 
10 11 Mq) star-forming galaxies since z « 1. Taken to- 
gether, these results indicate that quenching — the rapid 
cessation of star formation in galaxies — is an important 
driver of galaxy evolution at z < 1. Moreover, whatever 
mechanism or mechanisms are responsible for quenching 
must affect galaxies spanning a wide range of stellar mass 
at these redshifts. 

We can use the results presented in Section 5 to quan- 
tify both the quenching rate — the frequency with which 
star-forming galaxies are being transformed into quies- 
cent galaxies — and the stellar mass growth of the qui- 
escent galaxy population due to the addition of newly 
quenched galaxies. The evolution with redshift of both 
these quantities should place important constraints on 
the broad range of proposed quenching processes, and 
their implementation into theoretical galaxy formation 
models. 

In the subsequent analysis we make the simplifying but 
well-motivated assumption (see Section 6.1) that we can 
neglect the effects of mergers. More specifically, we im- 
plicitly assume that if mergers are taking place, they do 
not cause a significant number of galaxies to shift from 
one 0.5 dex wide stellar mass bin to another (e.g., from 
the 10 10 - 10 10 ' 5 M Q to the 10 10 ' 5 - 10 11 M Q stellar 
mass bin). Consequently, although we do not explicitly 
include their potential contribution, mergers could still 
be a viable means of quenching star formation in some 
galaxies. In addition, wc implicitly ignore in situ star 
formation within the quiescent galaxy population, which 
by design is negligible (see Figure 1). 

With the preceding caveats in mind, wc define the frac- 
tional quenching rate as -F qU cnch = (dnQ/dt)/rigF, or the 
fraction of star-forming galaxies that must be quenched 
per Gyr as a function of stellar mass and redshift in order 
to account for the measured evolution of the quiescent- 
galaxy population with decreasing redshift. Using equa- 
tion (9), we obtain 

F — «fe) ,1 + 2,,Q " ,SF "(5)"' (13) 
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Fig. 16. — Fractional quenching rate -F^e^h vs. redshift in four intervals of stellar mass between 10 95 - 10 115 Mq. F quonch is the 
proportion of star-forming galaxies in a given stellar mass interval that must be quenched per Gyr in order to match the measured evolution 
in the number density of quiescent galaxies. Although we implicitly ignore the effects of mergers in calculating F^e,^, that does not 
necessarily mean that mergers are not responsible for quenching star formation in some galaxies (see the discussion in Section 6.2). Although 
the uncertainties are large, we find that F quonc j 1 is typically a factor of ~ 2 — 3 higher for 10 9 ' 5 — 10 10 5 Mq galaxies compared to more 
massive galaxies > 10 11 Mq. Moreover, the fractional quenching rate in lower-mass systems appears to be increasing toward the current 
epoch. 



where the Q and SF subscripts refer to quiescent and 
star-forming galaxies, respectively, and the relevant co- 
efficients and uncertainties are listed in Table 7. The 
(dt/dz) term is the derivative of the age- redshift func- 
tion, t(z), where t is the age of the Universe (Hogg 1999). 
For convenience, we note that in our adopted cosmology 
(dt/dz) can be approximated from < z < 1 to an ac- 
curacy of better than 0.5% by a third-order polynomial 
of the form (dt/dz) -13.8835 + 19.3598z - 13.621z 2 + 
4.2141z 3 Gyr" 1 . 

In Figure 16 we plot .Fquench versus redshift in four 
bins of stellar mass between 10 9 5 — 10 115 Mq. This 
figure reveals several interesting results. First, we find 
that -French is relatively small at all redshifts and stellar 
masses; it varies between sa 5% — 10% Gyr -1 , except in 
the lowest stellar mass bin at z < 0.2, where it rises to 
« 12% - 18% Gyr" 1 . For comparison, Blanton (2006) 
estimate that roughly 25% of blue star-forming galax- 
ies must be quenched since z sa 1 in order to match 
the measured buildup of the optical luminosity func- 
tion for red-sequence galaxies. Second, although the 
uncertainties are large, our measurements suggest that 
-Fquench varies systematically with stellar mass. Among 
> 10 11 Mq galaxies, we find a typical quenching rate 
of « 5% Gyr -1 , whereas among < 10 10 " 5 Mq galax- 
ies the rate is around sa 7% — 15% Gyr -1 , a factor of 
^2 — 3 higher. Finally, Figure 16 suggests that F quenc h 
increases toward low redshift among lower-mass galax- 



ies. Specifically, we find a weak trend of an increasing 
-Fquonch toward z sa among 10 10 — 10 10 5 Mq galax- 
ies, and a much more significant and rapid rise toward 
low redshift among 10 95 — 10 10 Mq galaxies. Although 
the uncertainties are large, these results suggest that star 
formation quenching is more prevalent among lower-mass 
galaxies, and that the fraction of low-mass star-forming 
galaxies that are being quenched is increasing toward the 
current epoch. 

We can extend this analysis one step further and cal- 
culate Psf-s-q = d/9Q/dt, the rate of stellar mass transfer 
from the population of star-forming to quiescent galaxies, 
using the measured stellar mass growth of the quiescent- 
galaxy population. Once again, we implicitly assume 
that the stellar mass growth of the quiescent popula- 
tion is entirely due to the addition of newly quenched 
(i.e., previously star- forming) galaxies, and that dissipa- 
tionless mergers between two quiescent galaxies does not 
cause a significant number of galaxies to change their 
stellar mass by more than a factor of three. With these 
caveats in mind, we use equation (10) to write /?sf^q in 
units of Mq yr -1 Mpc -3 as 

Psf^q = Po,Q /8q (1 + zfQT 1 (J^j , (14) 

where the relevant coefficients are given in Table 7. 
In Figure 17 we plot psF-s-Q versus stellar mass in four 
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Fig. 17. — Rate at which stellar mass is being transferred from the star-forming to the quiescent galaxy population, psp_>Q, vs. stellar 
mass in four intervals of redshift between 2 = — 0.8. The individual symbols with error bars (offset slightly horizontally for clarity), 
correspond to different redshift bins, and the grey shaded region shows the mean trend and la confidence region over all redshifts. We find 
that pSF— »Q is largely independent of redshift, and peaks around ~ 10 10 8 Mq. Consequently, although a larger proportion of lower-mass 
galaxies are being quenched (see Figure 16), the bulk of the stellar mass buildup of the quiescent galaxy population is occuring near the 
'knee' of the SMF. 



equal-sized bins of redshift from z = — 0.8. The indi- 
vidual symbols with error bars, which have been offset 
slightly in the horizontal direction for clarity, correspond 
to different redshift bins, and the grey shaded region re- 
flects the broad trend we deduce. We find that psf->q 
depends weakly on redshift, but has a very strong stel- 
lar mass dependence. Quantitatively, /0sf->q rises from 
m 1.5 x 1(T 3 Mq yr" 1 Mpc -3 around ~ 10 98 Mq to 
a peak value of « 3 x 10 -3 Mq yr _1 Mpc -3 around 
~ 10 10 - 8 Mq. Above - 10 11 M , psf^q declines 
sharply to a mean value of«lx 10~ 3 Mq yr _1 Mpc -3 . 

These results reveal that although fractionally more 
low- and intermediate-mass (~ 10 —10 105 Mq) galax- 
ies are being quenched (see Figure 16), the bulk of the 
stellar-mass buildup within the quiescent galaxy popula- 
tion occurs around ~ 10 108 Mq, near the 'knee' of the 
SMF. As emphasized by Bell et al. (2007), this stellar 
mass scale appears to be important because both the stel- 
lar mass- weighted star-formation rate density in galaxies, 
and the stellar mass growth of the quiescent galaxy pop- 
ulation peak around this stellar mass at z < 1; however, 
the underlying physical cause of this coincidence remains 
unknown. 

7. SUMMARY 

Wc have measured the evolution of the SMFs of quies- 
cent and star-forming galaxies from z = — 1 using two 
large, statistically complete, spectroscopic samples. At 
low redshift we use a sample of ~ 170, 000 SDSS galaxies 



with GALEX, 2MASS, and WISE photometry, and at 
intermediate redshift we use a sample of ~ 40, 000 galax- 
ies brighter than i ss 23 drawn from PRIMUS with deep 
GALEX and IRAC imaging. Our PRIMUS sample is 
notable for its depth, sample size, and area, which totals 
« 5.5 deg 2 over five widely-separated fields, while our 
SDSS- GALEX sample comprises one of the largest sta- 
tistical samples of local galaxies with SDSS and GALEX 
photometry ever assembled. 

The exceptional multi-wavelength coverage of both 
datasets provides deep UV to mid-infrared imaging 
over the entire spectroscopic survey area, allowing us 
to robustly estimate SFRs and stellar masses using a 
new Bayesian SED-modeling code (iSEDfit; see Ap- 
pendix A). We use these measurements to separate the 
galaxy population into quiescent and star-forming based 
on their position in the SFR-stcllar mass diagram, and to 
measure the evolution of the SMF over a large dynamic 
range in stellar mass and redshift with relatively small 
sample variance and Poisson uncertainties. In addition, 
we carefully assess the effect of systematic errors in our 
stellar mass and SFR estimates, and find that the evolu- 
tionary trends we measure are broadly insensitive to the 
exact choice of priors and population synthesis models. 

Our principal quantitative results are as follows: 

1. We find for the global galaxy population that 
the SMF has evolved relatively little since z — 
1, although we do find evidence for differential 
evolution — mass assembly downsizing. We mea- 
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sure a 31% ± 9% increase in the integrated number 
density of ~ 10 10 Mq galaxies since z ps 0.6, and 
a —3% ± 6% change in the number density of all 
~ 10 11 A4q galaxies since z ps 0.9. Most massive 
galaxies, therefore, appear to be largely in place 
by z = 1, while lower- mass galaxies continue to 
assemble toward the present epoch. 

2. The relatively subtle changes in the SMF of the 
global population, however, mask much more dra- 
matic evolution in the SMFs of star-forming and 
quiescent galaxies. Within the star-forming pop- 
ulation the most rapid evolution occurs among 
massive galaxies, whereas the low-mass end of the 
star-forming galaxy SMF does not change signifi- 
cantly. We find that the comoving number density 
of 10 — 10 11 A4® star-forming galaxies changes by 
—8% ± 10% between z ss 0.8 and z ~ 0.0, whereas 
the space-density of massive (10 11 — 10 115 M.q) 
star-forming galaxies declines by 54% ± 7% since 
z w 1. 

3. Meanwhile, among quiescent galaxies the most sig- 
nificant evolutionary changes occur among low- 
and intermediate-mass (10 9 5 — 10 10,5 A4@) galax- 
ies, whereas most massive (> 10 11 A4®) quiescent 
galaxies are largely in place from 2 = — 1. Quan- 
titatively, we find a factor of ~ 2 — 3 increase in 
the number-density of 10 95 — 10 10 A4q galaxies 
since z w 0.5, and a marginally significant increase, 
22%± 12%, in the space-density of > 10 11 M Q qui- 
escent galaxies since z ~ I. 

We use these measurements to place new constraints 
on the growth of galaxies by mergers, and to quantify 
the buildup of the quiescent galaxy population due to 
star formation quenching as a function of redshift and 
stellar mass: 

1. Using a simple model to account for the expected 
growth of galaxies due to star formation, we find 
that mergers do not appear to be a dominant chan- 
nel for the stellar mass buildup of galaxies at z < 1 , 
even among massive (> 10 11 A4q) systems. Quan- 
titatively, we find that mergers are responsible for 
a net change in the number density of 10 9 ' 5 — 
10 10 - 5 M Q galaxies of-13%±9% Gyr" 1 from z 
to z « 0.5, and a change of -14% ± 18% Gyr -1 in 
the number density of 10 10 5 — 10 115 M Q galax- 
ies since z « 0.8. These results do not imply that 
mergers are not occurring, only that they do not 
have a large net effect on the SMF over the range 
of stellar masses and redshifts probed by PRIMUS. 

2. Our results also imply that the rate at which star 
formation is quenched in galaxies depends both 
on stellar mass and redshift, with a peak around 
~ 10 10 ' 8 A4q, and an increase in the quenching 
rate at lower redshift for lower mass galaxies. In 
particular, we find that the quenching rate for mas- 
sive galaxies with > 10 11 M.q is consistently low 
(f« 5% Gyr" 1 ) at all redshifts to z ps 1. 

To fully characterize the build up of stellar mass for all 
galaxies to z = 1, additional measurements at the lower 



mass end of the SMF are needed, which requires both 
deeper multi-wavelength imaging and spectroscopy. In 
addition, while our results at intermediate redshift use 
a large sample of ~ 40, 000 galaxies across five separate 
fields covering ps 5.5 deg 2 , they are still dominated by 
sample variance. This result argues that even more wide- 
field imaging and spectroscopy to at least i = 23 are 
needed to precisely measure the SMF from z = — 1. 
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APPENDIX 

A. iSEDf it SPECTRAL ENERGY DISTRIBUTION MODELING CODE 

A.l. Background 

iSEDf it 39 was developed in the idl programming language to be a fast and flexible tool to extract the physical 
properties of both nearby and high-redshift galaxies from their broadband UV, optical, and near-infrared SEDs. 
It builds on the Bayesian formalism pioneered by Kauffmann et al. (2003a) to model the optical spectral features 
of SDSS galaxies, and its subsequent extension to broadband photometry (e.g., Bundy et al. 2005; Salim et al. 2005, 
2007; da Cunha et al. 2008; Auger et al. 2009; McGee et al. 2011; Taylor et al. 2011). To date, stellar masses and other 
physical parameters derived using iSEDf it have been used to measure the evolution of the stellar mass-mctallicity 
relation (Moustakas et al. 2011), the stellar mass dependence of AGN accretion (Aird et al. 2012), the stellar mass 
and SFR surface densities of compact starbursts with high-velocity outflows (Diamond-Stanic et al. 2012), and the 
ages, SFRs, and stellar masses of galaxies at z = 6 — 10 (Zitrin et al. 2012; Zheng et al. 2012), with many subsequent 
applications forthcoming. 

Like all SED-modcling codes designed to extract the physical properties of galaxies 40 , iSEDfit relies on stellar 
population synthesis (SPS) models to provide as input S(X,t, Z), the spectral evolution of a simple stellar population 
(SSP) of a given stellar metallicity, Z. An SSP is an idealized stellar population formed in an instantaneous burst of 
star formation which evolves passively thereafter with time t. The basic ingredients of an SSP are (e.g., Tinslcy 1968; 
Bruzual 1983): (1) stellar evolution calculations for stars spanning the full range of initial mass; (2) a stellar library 
that provides the emergent spectrum of a star at each position in the Hertzsprung-Russell (HR) diagram; and (3) an 
assumed IMF, which specifies the relative number of stars of a given stellar mass. 

Unfortunately, even the relatively "simple" goal of modeling SSPs is limited by uncertainties in calculating particular 
phases of stellar evolution, inadequecies in the stellar libraries (e.g., non-solar abundance ratios, spectral coverage and 
resolution, etc.), and other simplifying assumptions (see Conroy et al. 2009 and Conroy & Gunn 2010 for recent in- 
depth discussions) . For example, among the least well understood phases of stellar evolution are the thermally pulsating 
asymptotic giant branch (TP-AGB) stars, blue stragglers (BS), and horizontal branch (HB) stars, all of which are 
relatively luminous and can therefore affect the integrated spectrum of the stellar population (e.g., Maraston 2005; 
Melbourne et al. 2012). SSP calculations also implicitly assume a well-sampled (i.e., fully populated), IMF, which may 
not always be satisfied (e.g., Fumagalli et al. 2011). 

Differences in how these issues are addressed (or ignored) for a given IMF can lead to significant systematic dis- 
crepancies among SSPs derived using different SPS models. In principal, the uncertainties affecting SSPs should be 
incorporated into the SED modeling in order to obtain reliable parameter estimates and realistic confidence intervals 
(e.g., Conroy et al. 2010). In practice, however, this procedure is both cumbersome and computationally challenging. 
Instead, iSEDf it adopts the simplified approach of allowing the user to select from among many different SPS models, 
thereby allowing the effects of choosing one set of SSPs over another to be systematically investigated. 

A. 2. From Simple Stellar Populations to Spectral Energy Distributions of "Galaxies" 

Given S(X, t, Z), iSEDf it computes the integrated SED of a "galaxy" (a composite stellar population) as a function 
of time t using the following convolution integral: 

C{\,t,Z)= [ 4it-t')S[X,t',Z(t-t')}lQ-° AA ^ t ">dt', (Al) 
Jo 

39 https : //github. com/moustakas/iSEDf it 40 ..^ j^-^- j ™r i i, ± ^ /omi\ 

v 6 u Sec http://www.sedfittmg.org and Walcher et al. (2011) 

for extensive references. 



Evolution of the Stellar Mass Function 



27 



\ 



o 



iafl 
o 



0.4 
0.2 
0.0 
-0.2 
-0.4 

0.4 
0.2 
0.0 
-0.2 



K— correct - 
: N= 165473 : 


'- MPA/JHU-DR4 1 
: N=123623 : 




- \ - 


'- MPA/JHU-DR7 : 

N=161527 : 


Salim + 07 - 

: N=34056 : 







10 11 

log (M / M ) iSEDfit 



12 9 10 11 

log (M / M ) iSEDfit 



12 



Fig. Al. — Comparison of stellar masses derived using iSEDfit (see Section 4.1) for our SDSS-GALEX sample vs. stellar masses derived 
using the four independent techniques described in Appendix A. 4. The solid line in each panel indicates the median residual trend, and the 
dashed lines show the interquartile range of the residuals. The dotted horizontal line shows for reference the one-to-one relation. Overall, 
we find very good agreement between the stellar masses derived using iSEDf it and these various other methods. 

where ip(t) is the star-formation history (SFH), and A(X) is the wavelength-dependent attenuation, which in general 
may depend on time (e.g., Chariot & Fall 2000). The current version of iSEDf it only handles mono-metallic stellar 
populations (i.e., it does not treat the chemical evolution of the system self-consistently) in which case equation (Al) 
reduces to 



C(X,t,Z) = / ^(t- 
Jo 



t')S(X,t',Z) io-° AA ^dt'. 



(A2) 



In order to solve this integral, iSEDfit requires several additional inputs (implicit prior assumptions) from the 
user. First, an extinction or attenuation 41 curve, fc(A) = A(X)/E(B — V), must be chosen among several different 
possibilities, where A(X) is the total wavelength-dependent attenuation and E(B — V) is the color excess (Calzetti 
2001). The currently supported possibilities are the Calzetti et al. (2000) starburst galaxy attenuation curve, the 
Chariot & Fall (2000) attenuation law, the O'Donnell (1994) Milky Way extinction curve, the extinction curve of 
the Small Magellanic Cloud (Gordon et al. 2003), or no attenuation. Note that among these only the Chariot & Fall 
(2000) attenuation curve is time-dependent. 

Next, a parametric form for the SFH must be specified. In principal, ip(t) could be arbitrarily complex, or even 
non-parametric (e.g., Cid Fernandcs et al. 2005; Tojeiro et al. 2007; Panter et al. 2008). In practice, however, it is 
challenging if not impossible to constrain the detailed SFHs of individual galaxies from broadband photometry alone 
(i.e., without high- resolution spectroscopy, although see Dye 2008). On the other hand, the integrated spectra of 
many star- forming galaxies may be poorly fitted by simple (e.g., exponentially declining) SFHs. Therefore, iSEDfit 
optionally allows stochastic bursts to be superposed on a backbone of smooth SFHs (e.g., Kauffmann et al. 2003a). For 
the underlying smooth component, the user can choose either exponentially declining SFHs (so-called simple r-models; 
Sandage 1986): 



M 



tot e _t/ T . 



(A3) 



or "delayed" r-modcls: 



41 Recall that attenuation includes the effects of both absorp- 
tion and scattering, whereas extinction describes the absorption of 



light by a homogenous foreground dust screen (Witt et al. 1992; 
Witt & Gordon 2000). 
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M 



Mt) = ^te-^\ (A4) 

where the subscript "s" indicates that these are "smooth" SFHs, t is the age of the stellar population (the time since 
the onset of star formation), r is the characteristic time for star formation, and the normalization is defined to be 
■Mtot = 1 -Mq- The delayed r-modcls are advantageous because they allow for both exponentially declining (t/r ^> 1) 
and linearly rising (t/r <C 1) SFHs to be explored, the latter of which are needed to accurately reproduce the colors 
of high-redshift (z > 2) galaxies (e.g., Maraston et al. 2010; Papovich et al. 2011; Behroozi et al. 2012). 

We characterize each burst by three independent parameters: the time the burst begins, tf>, its duration, Atj,, and 
the burst fraction, the relative strength of the burst, J-f,. The SFH of each burst, ipb(t), is a Gaussian function given 
by 

V , b (i) = ^ e -('-M 2 /2A^ (A5) 

where Ab is the burst amplitude. Defining Th to be the mass formed in the burst divided by the total mass formed by 
the underlying r model until the peak of burst, we obtain 

T - Mb Atb A ( AP.} 

J ^W = AMT^j A ' (A6) 

where A4 T (tb) = J^ b tp T (t)dt, and Adtot = 1 M.@. Equation (A6) assumes a simple r-model, but a similar expression 
can be derived for the delayed r-model. The final composite SFH is given by 

m = Mt)+ E ^(t), (A7) 

where A^urst is the total number of bursts experienced by each model galaxy. The number of bursts is determined by 
specifying the cumulative probability Pburst for a burst to occur within a APburst time interval. Finally, we note that 
the current version of iSEDf it additionally allows the final burst to be truncated exponentially with a characteristic 
time r trunc , thereby allowing the SFHs and physical properties of post-starburst galaxies to be investigated (e.g., 
Tremonti et al. 2007). 

A. 3. Extracting the Physical Properties of Galaxies from Broadband Photometry 

Based on the large number of free parameters needed to model the integrated SEDs of galaxies, it would be far 
too computationally expensive to explore all possible parameter combinations (e.g., on a uniform grid). Moreover, 
traditional best-fitting (maximum likelihood) techniques arc limited because they only account for photometric uncer- 
tainties, but not physical degeneracies among different models (parameter combinations). Alternatively, Markov Chain 
Monte Carlo (MCMC) algorithms may be more suitable for exploring the multi-dimensional parameter space (e.g., 
Acquaviva et al. 2011). However, MCMC methods are typically too slow to enable tens, hundreds, or even millions of 
galaxies to be fitted in a timely manner with multiple independent prior parameter combinations or SPS models. 

Given the limitations of these other techniques, iSEDf it extracts the physical parameters of interest using a simplified 
Bayesian approach (Walcher et al. 2011, and references therein). First, the model parameters are drawn from a user- 
specified prior probability distribution using a Monte Carlo technique. Next, given the broadband fluxes F± of a galaxy 
at redshift z in i = 1, N filters, and the corresponding cr^ uncertainties, iSEDfit uses Bayes' theorem to compute the 
posterior probability distribution function (PDF) 

p(Q\F i ,z)=p(Q)xp(F i ,z\Q), (A8) 

where Q represents the set of model parameters (stellar mass, age, mctallicity, etc.). Here, p(Q) is the prior probability 
of the model parameters, p(Fi,z\Q) is the likelihood C cx exp[— x 2 (Fi, z|Q)/2] of the data given the model, and x 2 is 
the usual goodness-of-fit statistic appropriate for normally distributed photometric uncertainties, given by 

1 = 1 Gl 

where A is a normalization factor, and the C^(Q,z) are the broadband fluxes of each model SED given the redshift 
and parameter combination Q. Once % 2 has been computed for every model, the marginalized posterior PDF of the 
parameter of interest, for example p(A4) for the stellar mass, can be derived by randomly drawing each parameter 
value with probability given by equation (A8), thereby effectively integrating (using histogram binning) over the 
other "nuisance" parameters. Although iSEDf it is capable of reconstructing the full posterior distributions in post- 
processing after the computationally intensive fitting has been completed, the code by default provides the median 
of the posterior PDF as the best estimate of each parameter, and estimates the uncertainty in each parameter as 
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Fig. Bl. — Effect on the SDSS-GALEX SMF of (left) different prior parameter combinations and (right) different population synthesis 
models (see the text in Appendix B for acronym definitions and details). We find that assuming that galaxies do not experience stochastic 
bursts ("no bursts") has a significant effect on the massive end (A4 > 10 11 Mq) of the SMF, whereas for all other combinations of priors 
and population synthesis models the effects are relatively small (< O.f dex). 

1/4 of the 2.3 — 97.7 percentile range of the posterior distribution, which would be equivalent to la for a Gaussian 
distribution. 

A. 4. Consistency Checks 

In Figure Al we verify that the stellar masses we derive using iSEDfit for our SDSS-GALEX sample based on 
our fiducial prior parameters (see Section 4.1 and Appendix B) are reasonable, by comparing them against the stellar 
masses of the same objects derived using four other independent techniques. In each panel the solid line indicates the 
median residual trend, and the dashed lines show the interquartile range of the residuals. 

The upper-left panel compares our stellar masses against the stellar masses derived using K-correct 
(Blanton & Roweis 2007), fitted to the same 12-band UV to mid-infrared photometry as iSEDf it (see Section 2.4). We 
find a weak residual trend with stellar mass, in the sense that K-correct yields slightly higher (lower) stellar masses for 
lower-mass (higher-mass) galaxies, which for the massive galaxies at least is similar to the residual trend reported by 
Bernardi et al. (2010). The upper-right and lower- left panels compares our masses against the stellar masses publicly 
released by the MPA/JHU team based on the SDSS DR4 42 and DR7 43 data releases. The MPA/JHU-DR4 masses 
are based on fitting the US a and Z? n (4000) optical spectral indices (Kauffmann ct al. 2003a), and the DR7 stellar 
masses are derived using a similar technique as iSEDf it, but only fitting to the SDSS ugriz photometry. In both 
cases the agreement between the various independent determinations is outstanding. Finally, in the lower-right panel 
we compare our mass estimates against the stellar masses derived by Salim et al. (2007) using a similar SED-modcling 
technique as iSEDf it, but just fitting to the GALEX plus ugriz photometry. Once again, the agreement between the 
two stellar mass estimates is very good. 

B. EFFECT OF PRIOR PARAMETER CHOICES AND POPULATION SYNTHESIS MODELS ON 

OUR RESULTS 

In this appendix we examine the effect of varying the SPS models and prior parameters we use to derive stellar masses 
and SFRs on our results (see also Section 4.1). We consider four distinct SPS models: FSPS (v2.3; Conroy et al. 2009; 
Conroy & Gunn 2010); Bruzual & Chariot (2003, hereafter BC03); the SPS models of Maraston (2005, hereafter 
Ma05) 44 ; and pegase 45 (Fioc & Rocca-Volmerange 1997, 1999; Le Borgne et al. 2004). For the FSPS and BC03 
models we adopt the Chabrier (2003) IMF from 0.1 - 100 Mq, and for Ma05 and pegase we use the Kroupa (2001) 
IMF from 0.1 — 100 M.q. We neglect the ~ 0.03 dex systematic difference between the two IMFs. Each of these 



42 
43 



http : //www . mpa-gar ching . mpg . de/SDSS/DR4 44 



http://www-astro. physics, ox. ac.uk/~mar aston/Claudia's_Stellar_Populatic 
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Fig. B2. — Relative number density vs. redshift for (top) quiescent and (bottom) star-forming galaxies in four bins of stellar mass between 
10 9 ' 5 - 10 11 ' 5 M Q based on different SPS models and prior assumptions (see text). Each relation has been normalized by the number 
density at z = 0.1 so that the relative evolutionary trends can be compared. The grey shaded region reflects the la confidence region of 
the mean fitted relation based on our fiducial SPS models and prior assumptions (see Figure 12). We find that the individual evolutionary 
trends are generally consistent with one another at the ±1<t level except for the stellar masses and SFRs derived using the Calzetti et al. 
(2000) dust law, which exhibits a shallower (steeper) decline (increase) in the number of massive star-forming (quiescent) galaxies. Overall, 
however, we conclude that our results are broadly insensitive to our choice of SPS models and priors. 

models relies on a different combination of stellar evolution calculations and stellar libraries and therefore differ in 
their predictions of the integrated spectra of galaxies. 

For completeness, we briefly summarize the salient features of each of these SPS models. Our fiducial stellar masses 
and SFRs are based on the empirically calibrated version of the FSPS models described by Conroy & Gunn (2010). 
We couple these models to the Padova 46 stellar evolutionary isochrones (Girardi et al. 2000; Marigo & Girardi 2007; 
Marigo et al. 2008), which have been supplemented with the Vassiliadis & Wood (1994) models for post-AGB stars. 
Integrated spectra are generated using the low-resolution, semi-empirical BaSeL 3.1 library (Lejeune et al. 1997, 1998; 
Westera et al. 2002), which extends from the UV to the rest- frame near- infrared, except for the TP-AGB stars, 
for which the empirical spectra of Lancon & Mouhcine (2002) are used over the full wavelength range. The BC03 
models we use are based on the Padova (1994) isochrones (Alongi et al. 1993; Bressan et al. 1993; Fagotto et al. 1994), 
supplemented with the Vassiliadis & Wood (1993) and Vassiliadis & Wood (1994) models for TP-AGB and post-AGB 
stars, respectively. Integrated spectra are synthesized using the empirical STELIB stellar library (Le Borgne et al. 
2003) in the optical (3200 — 9500 A) and extended into the UV and near-IR at lower resolution using the BaSeL 3.1 
library. Finally, the Ma05 models are based on the stellar tracks and isochrones through the main-sequence turnoff 
published by Schaller et al. (1992) and Cassisi et al. (1997, 2000). The fuel consumption theorem (Maraston 1998, 
2005) is used to calculate the post-main sequence phases of stellar evolution, including the TP-AGB phase. These 
evolutionary calculations are then tied to the Langon & Mouhcine (2002) empirical spectra for TP-AGB stars and to 
the BaSeL 3.1 spectral library for other stellar populations, resulting in a set of integrated spectra from the UV to 
the near-infrared. In our calculations we adopt the version of the models computed using the "red horizontal branch" 
morphology (see Maraston 2005 for details). Finally, pegase utilizes the early-1990s version of the Padova stellar 
isochrones and couples those to the BaSeL stellar library. 

Our fiducial prior parameters were briefly described in Section 4.1, but here we provide more details (sec Appendix A 
for additional salient details and parameter definitions) . We assume exponentially declining SFHs with Gaussian bursts 
of varying onset, strength, and duration randomly superposed. Following Kauffmann et al. (2003a) and Salim et al. 
(2007), we draw t^ 1 from a uniform distribution in the range 0.01 — 10 Gyr -1 and allow bursts to occur with a 
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cumulative probability Pburst = 0.5 every APburst = 2 Gyr. We draw Atj, from a logarithmic distribution in the range 
30 — 300 Myr (i.e., shorter-duration bursts are preferred) and Fb from a logarithmic distribution spanning 0.03 — 4 
(Salim ct al. 2007; Wild et al. 2009). We allow the age t (time for the onset of star formation) of each model to 
range with equal probability between 0.1 — 13 Gyr, although we disallow ages older than the age of the Universe at 
the redshift of each galaxy. We assume a uniform prior on stellar metallicity Z in the range 0.004 — 0.03 (roughly 
20% — 150% times the solar metallicity; Asplund et al. 2009). Because the SSP models are generally only available 
for a small number of tabulated values of Z, we linearly interpolate between these values to obtain an SSP with 
an arbitrary metallicity. Finally, we adopt the time-dependent attenuation curve of Chariot & Fall (2000), in which 
stellar populations older than 10 Myr are attenuated by a factor /i times less than younger stellar populations. We 
draw /i from an order four Gamma distribution that ranges from zero to unity centered on a typical value (fi) =0.3 
(Chariot & Longhetti 2001; Wild et al. 2011) and the U-band optical depth from an order two Gamma distribution 
that peaks around Ay sa 1.2 mag, with a tail to Ay ss 6 mag. 

We consider the effect of varying a small number of these priors and SPS models on our results. Specifically, we 
consider stellar masses and SFRs derived assuming: (1) the Calzetti et al. (2000) starburst galaxy attenuation curve 
("Calzetti"); (2) that galaxies do not experience stochastic bursts of star formation ("no bursts"); and (3) fixed solar 
metallicity ("solar metallicity"). Although these parameter combinations are not exhaustive, they have been chosen to 
reasonably span the range of priors commonly adopted in other studies of the SMF (see, e.g., Marchesini et al. 2009). 

In Figure Bl we plot the SDSS- GALEX SMF derived using these different priors and SPS models. In the left panel 
we use our fiducial FSPS models and vary the prior assumptions, and in the right panel we use the same fiducial 
priors and vary the SPS models. Overall we find that these variations have a relatively small systematic effect on the 
SMF. The most significant differences result when we do not include the effects of bursts, which leads to typically 
higher stellar masses for massive galaxies. In Figure B2 we plot number density versus redshift for (top) quiescent and 
(bottom) star-forming galaxies in four bins of stellar mass between 10 9 5 — 10 115 M.q (see Section 5.2). For clarity 
we only show the slope of the line [in log-log space; see equation (9)] fitted to the mean number density of galaxies 
measured in each redshift bin, normalized by the number density at z = 0.1 so that the relative evolutionary trends can 
be compared. In each panel, the grey shaded region reflects the la confidence region of the mean fitted relation based 
on our fiducial SPS models and prior assumptions (see Figure 12 and Table 7). We find that the relative evolutionary 
trends we infer are generally within ±lcr of the trends inferred using these other SPS models and prior assumptions. 
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TABLE 1 

Sample Properties 





Selection 


Magnitude 
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Limit a 


(deg 2 ) 


N c 


CDFS d 
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23.0 


1.496 


8050 
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7290 


ELAIS-Sl d 
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23.2 


0.800 


4140 


XMM-SXDS 
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23.0 


0.646 


6403 


XMM-CFHTLS 


i 


23.0 


1.700 


14547 


PRIMUS" 






5.499 


40430 


SDSS- GALEX 


r 


17.6 


2504 


169727 



a Corrected for foreground Galactic extinction. 
b Angular area surveyed. 

c Number of galaxies satisfying the selection criteria given in 
§3.1. 

d Due to differences in the PRIMUS experimental design, 
our CDFS and ELAIS-S1 samples arc also flux-limited in the 
3.6 /im IRAC band between 17 < [3.6] < 21. 
° Combination of all five fields. 



TABLE 2 
Stellar Mass Completeness Limits' 1 





COSMOS 


XMM-SXDS XMM-CFHTLS 


CDFS 


ELAIS-S1 


Redshift Range 




log (Mi 


imAM©) 






All 


0.20 - 0.30 


8.73 


8.86 


8.95 


9.62 


9.70 


0.30 - 0.40 


9.14 


9.23 


9.23 


9.87 


9.99 


0.40 - 0.50 


9.51 


9.58 


9.51 


10.10 


10.26 


0.50 - 0.65 


9.92 


9.97 


9.87 


10.37 


10.56 


0.65 - 0.80 


10.33 


10.38 


10.31 


10.65 


10.87 


0.80 - 1.00 
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Star-Forming 


0.20 - 0.30 
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9.58 
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Quiescent 
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a Stellar mass completeness limits among all, quiescent, and star-forming galaxies as a 
function of redshift. Above these limits our sample includes more than 95% of all types 
of galaxies, accounting for the flux limit in each field and mass-to-light ratio variations. 
For comparison, in our SDSS- GALEX sample the completeness limit is 10 9 Mq for 
all three samples. 
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SDSS-GALEX Stellar Mass Function 
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TABLE 4 

PRIMUS Stellar Mass Function for. All, Quiescent, and 
Star-Forming Galaxies 
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Note. - PRIMUS stellar mass function at z = 0.2 — 1 for all, star-forming, and quiescent 
galaxies. The units of $ are 10~ 4 Mpc~ 3 dex -1 , a cv is the estimated la uncertainty in $ due to 
sample variance, and iV is the number of galaxies in each stellar mass bin. 
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TABLE 5 

Cumulative Number and Stellar Mass Density of All Galaxies 



<*) 


log (n) a 


log(p) b 


log (n) 


log (p) 


log (n.) 


log (p) 


log (n) 


log (p) 




log (M/M 


) >9.5 


log (M/M 


©) > 10 


log (M/M 


©) > 10.5 


log (M/M 


0) > 11 


0.100 


-2.09 ±0.05 


8.35 ±0.05 


-2.32 ±0.05 


8.31 ±0.05 


-2.68 ±0.05 


8.18 ±0.05 


-3.41 ±0.05 


7.78 ± 0.04 


0.250 


-2.12 ±0.05 


8.32 ± 0.09 


-2.36 ±0.06 


8.28 ±0.10 


-2.70 ±0.09 


8.16 ±0.13 


-3.41 ±0.29 


7.78 ±0.29 


0.350 


-2.11 ±0.06 


8.35 ± 0.06 


-2.30 ±0.05 


8.32 ±0.06 


-2.65 ±0.07 


8.20 ± 0.07 


-3.39 ± 0.10 


7.77 ±0.11 


0.450 


-2.16 ±0.09 


8.31 ±0.08 


-2.35 ±0.07 


8.27 ± 0.08 


-2.69 ±0.08 


8.16 ±0.09 


-3.47 ± 0.14 


7.70 ±0.16 


0.575 


> -2.22 


> 8.30 


-2.39 ±0.04 


8.28 ±0.05 


-2.69 ±0.04 


8.18 ±0.05 


-3.42 ±0.08 


7.78 ±0.10 


0.725 






> -2.37 


> 8.35 


-2.61 ±0.05 


8.27 ± 0.05 


-3.28 ±0.07 


7.90 ±0.08 


0.900 










> -2.79 


> 8.15 


-3.35 ± 0.07 


7.85 ± 0.07 



a Number density in h? Mpc 3 . 

b Stellar mass density in /170 Mq Mpc 



TABLE 6 

Number and Stellar Mass Density of All, Quiescent, and 
Star-Forming Galaxies 







log (n) 






log (p) 




(*) 




(h% Mpc" 3 ) 




(h 


70 M & Mpc" 










9.5 < log (M/M, 


s) < 10 








All 


Q" 




All 




SI' 


0.100 


-2.48 ± 0.05 


-2.94 ±0.07 


-2.66 ± 0.04 


7.28 ± 0.05 


6.83 ± 0.07 


7.10 ± 0.04 


0.250 


-2.49 ± 0.03 


-3.15 ± 0.04 


-2.60 ± 0.03 


7.27 ±0.03 


6.63 ± 0.05 


7.16 ± 0.03 


0.350 


-2.54 ±0.09 


-3.21 ±0.16 


-2.65 ±0.07 


7.24 ± 0.08 


6.58 ±0.13 


7.13 ±0.06 


0.450 


-2.61 ±0.15 


> -3.49 


-2.67 ±0.14 


7.16 ±0.13 


> 6.33 


7.10 ±0.12 


0.575 


> -2.71 




> -2.74 


> 7.06 




> 7.02 








10 < [og(M/M,^ 


) < 10.5 








All 


Q 


SF 


All 


Q 


SF 


0.100 


-2.56 ±0.05 


-2.89 ±0.05 


-2.84 ±0.04 


7.71 ±0.05 


7.39 ±0.05 


7.42 ± 0.04 


0.250 


-2.62 ±0.03 


-2.99 ±0.05 


-2.85 ±0.04 


7.65 ± 0.03 


7.29 ± 0.04 


7.39 ± 0.04 


0.350 


-2.56 ±0.04 


-3.01 ±0.04 


-2.75 ±0.04 


7.71 ± 0.03 


7.28 ± 0.04 


7.51 ±0.04 


0.450 


-2.62 ±0.04 


-3.14 ±0.07 


-2.78 ±0.03 


7.65 ±0.04 


7.16 ±0.05 


7.47 ± 0.04 


0.575 


-2.68 ±0.05 


> -3.20 


-2.84 ±0.04 


7.60 ± 0.05 


> 7.12 


7.42 ± 0.04 


0.725 


> -2.74 




> -2.88 


> 7.56 




> 7.40 








10.5 < log (M/M 


0) < 11 








All 


Q 


SF 


All 


Q 


SF 


0.100 


-2.77 ±0.05 


-2.99 ±0.05 


-3.17 ±0.04 


7.97 ±0.05 


7.76 ± 0.05 


7.55 ±0.04 


0.250 


-2.80 ±0.04 


-3.06 ± 0.04 


-3.15 ±0.07 


7.93 ± 0.04 


7.70 ± 0.04 


7.55 ± 0.09 


0.350 


-2.74 ±0.06 


-3.02 ± 0.09 


-3.07 ±0.04 


7.99 ± 0.06 


7.73 ± 0.08 


7.66 ± 0.04 


0.450 


-2.76 ±0.06 


-3.11 ±0.06 


-3.03 ±0.06 


7.97 ±0.06 


7.64 ±0.06 


7.69 ± 0.06 


0.575 


-2.78 ±0.03 


-3.11 ±0.04 


-3.06 ±0.03 


7.96 ± 0.03 


7.65 ±0.04 


7.66 ± 0.03 


0.725 


-2.72 ±0.05 


-3.06 ±0.04 


-2.99 ±0.07 


8.03 ± 0.05 


7.70 ± 0.04 


7.75 ± 0.07 


0.900 


> -2.93 


> -3.23 


-3.23 ± 0.07 


> 7.84 


> 7.56 


7.52 ±0.08 








11 < log(M/M 


) < 11.5 








All 


Q 


SF 


All 


Q 


SF 


0.100 


-3.42 ±0.04 


-3.54 ±0.04 


-4.04 ± 0.04 


7.74 ± 0.04 


7.63 ± 0.04 


7.09 ± 0.04 


0.250 


-3.42 ±0.20 


-3.55 ±0.24 


-4.10 ±0.17 


7.75 ±0.21 


7.60 ±0.27 


7.03 ± 0.28 


0.350 


-3.40 ± 0.04 


-3.55 ±0.06 


-3.94 ±0.03 


7.76 ± 0.04 


7.62 ±0.07 


7.20 ± 0.05 


0.450 


-3.48 ±0.10 


-3.71 ±0.11 


-3.88 ±0.09 


7.68 ±0.12 


7.45 ± 0.12 


7.27 ±0.10 


0.575 


-3.43 ±0.04 


-3.62 ±0.04 


-3.90 ±0.07 


7.74 ±0.05 


7.56 ± 0.04 


7.25 ± 0.08 


0.725 


-3.29 ±0.05 


-3.51 ±0.05 


-3.69 ±0.05 


7.88 ± 0.05 


7.66 ± 0.06 


7.47 ± 0.05 


0.900 


-3.37 ±0.05 


-3.58 ±0.03 


-3.77 ±0.08 


7.81 ±0.05 


7.59 ± 0.04 


7.40 ± 0.07 



Q = quiescent; SF = star- forming. 



Evolution of the Stellar Mass Function 



TABLE 7 

Power-Law Fits to the Number and Stellar Mass Density Evolution 51 



Samplc a 


log (n ) b 






All 

Q 

SF 




9.5 < log (M/Mq) < 10 




-2.441 ± 0.03 
-2.802 ± 0.03 
-2.666 ±0.04 


-0.729 ± 0.328 7.319 ± 0.03 
-3.444 ± 0.441 6.957 ± 0.04 
0.317 ±0.445 7.091 ±0.03 


-0.663 ±0.294 
-3.150 ±0.514 
0.351 ±0.400 


All 

o 

SF 




10 < log{M/M e ) < 10.5 




-2.527 ±0.04 
_o S21 + 03 
-2.830 ±0.04 


-0.644 ± 0.318 7.735 ± 0.04 
0.090 ±0.326 7.419 ±0.04 


-0.576 ±0.313 
_ i cci -|_ n 32fi 

0.168 ±0.339 


All 

Q 

SF 




10.5 < log (M/M e ) < 11 




-2.766 ± 0.02 
-2.956 ±0.02 
-3.150 ±0.06 


-0.062 ±0.145 7.968 ±0.02 
-0.775 ±0.144 7.786 ±0.02 
0.245 ±0.390 7.559 ±0.05 


-0.037 ±0.130 
-0.699 ±0.132 
0.344 ± 0.347 


All 

Q 
SF 




11 < log(M/M Q ) < 11.5 




-3.433 ± 0.02 
-3.534 ±0.02 
-4.090 ±0.01 


0.148 ±0.124 7.730 ±0.02 
-0.285 ±0.168 7.641 ±0.02 
1.114 ±0.100 7.036 ±0.01 


0.189 ±0.123 
-0.297 ±0.174 
1.274 ±0.109 



Q = quiescent; SF = star- forming. 

Model given by n(z) — no(l + z) 1 with n in h? Mpc -3 . 
Model given by p(z) = po(l + zy with p in /170 Mq Mpc 



